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ABSTRACT 


In this thesis, methods for forming a communieations link between Wireless 
Sensor Networks (WSNs) by enabling eaeh WSN to act as a smart antenna are presented. 
Each WSN is simulated as a set of randomly placed sensor nodes within a planar area. 
The proposed method involves a searching WSN, a receiving WSN and a link budget for 
establishing the link. 

The searching WSN has the task of transmitting a search beam in order to find 
adjacent WSNs. Like a lighthouse this is done in a rotating beam style search using the 
sensor nodes as an aperiodic array. Results show that for a random array, we can achieve 
a specific beamwidth and gain as a function of the number of elements and area. We also 
demonstrate that for a given required gain level we can spatially thin the array without 
significant loss of gain or the effects of grating lobes. 

The receiving WSN uses a spread spectrum based space division multiple access 
(SDMA) receiver. This receiver is simulated to determine the direction of arrival from the 
searching WSN and to extract the location information from the searching WSN’s signal 
in additive white Gaussian noise. From the DOA and the location information within the 
arriving signal, the WSN has sufficient knowledge to respond to the query of the 
searching WSN and form the communications link. 
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EXECUTIVE SUMMARY 


A wireless sensor node is a low cost, small, battery powered electronic device 
designed to monitor or measure a physical phenomenon of the environment around it, 
such as temperature, pressure, humidity, soil pH, vibration, motion, light, sound, 
radiation, and chemical presence. Wireless sensors are generally composed of a 
microprocessor, the appropriate sensor(s), a transceiver and a power source. Once the 
sensor makes a detection, the information collected is relayed back to an end user for 
further processing or analysis. A collection of sensors in close proximity to one another 
may form a Wireless Sensor Network (WSN) by establishing communications with one 
another through some form of self-organization. There are many different applications for 
WSNs, both commercially and militarily. In this work we assume a military application 
where the sensor nodes are deployed in hostile territory, and therefore delivery of the 
sensors may be via an aerial vehicle. Delivery in this manner would generally lead to a 
random spatial distribution of the sensors, yet the sensor nodes are close enough to form a 
local wireless sensor network. The nodes of this WSN may act collectively as a group 
for such actions as distributed processing or forming a distributed antenna array. 

In this thesis, we assume the deployed WSNs must transmit their information to 
the end user without access to the delivering aerial vehicle, but rather terrestrially. In the 
solution proposed here, where the WSNs are deployed to multiple locations nearby to one 
another, they can form a terrestrial Over The Horizon (OTH) network by forming a 
multi-hop communication link using the individual WSNs. In order to form this link, the 
WSNs must also solve the additional task of locating one another. The WSNs act 
collectively to form a smart antenna array. By establishing communication links between 
adjacent WSNs, an extended network made up of what has now become nodes of a larger 
Over The Horizon (OTH) network is formed. Once the OTH network is connected, the 
individual WSNs become part of a larger network that can relay sensed data to an end 
user. 

The goal of this thesis is to develop and validate methods to establish OTH 
communication links using wireless sensor networks by enabling each to act as a smart 



antenna array. In this thesis, we approach establishing the communication link between a 
searching WSN and a receiving WSN as a four-step process: localization of the sensor 
nodes within each WSN to determine its geometry; the search by one WSN for adjacent 
WSNs by acting as an smart antenna; determination of the direction of arrival and range 
of the searching WSN by the receiving WSN; and then forming the multi-hop link using 
WSNs by optimizing the gain of the beam forming the link. 

In searching for the adjacent WSN, each of the sensors act in concert as a smart 
antenna array. To describe how we will form a beam using a group of randomly placed 
sensor nodes, we review uniform linear and planar arrays then the random planar array. 
Each WSN is modeled as a uniformly random distribution of sensors nodes. In order to 
link the WSNs, we will wish to control the beamwidth, gain and direction of the array 
beam. We show analytically that we can steer the beam in an arbitrary direction. We also 
demonstrated control over the beamwidth and gain by providing MATLAB results and 
compare the results to analytical solutions. To demonstrate the search technique, a 40 
element, 2.5X x 2.5X random planar array scanning the horizon is presented. 

The results show that by choosing the physical size, likely as sub-section of the 
WSN we can choose the beamwidth of the search array. Within this area, we 
demonstrate control over the gain by choosing the number of elements employed in the 
array. Alternately, for a given required gain level we show that we can thin the array (or 
reduce density) dramatically without significant loss of gain or the effects of grating 
lobes. It does not matter which particular sensor nodes within the WSN are used as long 
as, on average, they are randomly distributed within the physical area chosen. This is 
important as we may wish to distribute the burden of transmitting, and thus consuming 
energy among the nodes within the WSN. 

From the perspective of the receiving WSN, this array needs to know the direction 
of arrival (DOA) and range to the searching WSN in order to form an optimum link. A 
new approach to DOA determination using spread spectrum space division multiple 
access (SDMA) is introduced, demonstrated for a uniform linear array and then validated 
for a random planar array. We show that the method works better for Walsh-Hadamard 



spreading codes than for pseudorandom noise codes and that 64 chips is sufficient to 
accurately establish the DOA and character of a linear antenna array. Using a random 
planar array, we determined the DOA of an incident test signal, then simulated extracting 
a message data stream from this signal in a noisy channel with SNRs of 20 db and -5 dB, 
accurately extracting the message. 

A link budget analysis is performed for a simulation involving two WSNs using 
random planar arrays on a flat earth. This example uses commercially available sensor 
motes at two transmit power levels. Gain calculations are made for each power level, 
solving for the required number of transmit and receive elements in each array. In this 
example we use the Crossbow MICAz® sensor operating at 2.4 GHz. The path loss 
provided the greatest amount of uncertainty in the link budget. Of the models examined, 
the range of values for a 500 m signal at 2.4 GHz varied by approximately 15 dB. After 
determining the gain of the random planar array, we found that for a 500 meter separation 
we could use 17 elements on the receive side and 25 elements on the transmit side when 
all elements were operating at full power. When all were operating at minimum power, 
we found that using 45 elements to receive and 45 elements to transmit provided 
sufficient margin. 

In this thesis we have developed and validated methods to establish 
communication links between wireless sensor networks using smart antenna 
beamforming techniques, spread spectrum space division multiple access as a method for 
a receiving WSN to determine the direction of arrival and message sent from an adjacent 
searching WSN and the use of link analysis to optimize the number of sensors used to 
form the link. 
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I. INTRODUCTION 


A wireless sensor node is a low eost, small, battery powered eleetronic device 
designed to monitor or measure a physical phenomenon of the environment around it, 
such as temperature, pressure, humidity, soil pH, vibration, motion, light, sound, 
radiation, and chemical presence. Wireless sensors are generally composed of a 
microprocessor, the appropriate sensor(s), a transceiver and a power source. Once the 
sensor makes a detection, the information collected is relayed back to an end user for 
further processing or analysis. Figure la shows an example of the MICAz sensor mote by 
Crossbow® Technologies [1] and in Figure lb a MICA2DOT sensor by Crossbow . 



(a) MICAz (b) MICA2DOT Sensor Mote 

Figure 1. Examples of Sensor Motes from Crossbow (From [1]) 


A collection of sensors in close proximity to one another may form a Wireless 
Sensor Network (WSN) by establishing communications with one another through some 
form of self organization. This group may act collectively to increase sensing or 
processing power by performing some of the computations onboard the sensor node itself 
[ 2 ]. 
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Many different types of applieations for wireless sensor networks exist, in both 
industry and the military world. Because they are wireless and self-organize, wireless 
sensor networks have the potential to relay information without a large infrastructure cost 
or physical impact. Commercial applications of wireless sensor networks include 
industrial monitoring, building controls, security, traffic management, weather, wildlife 
tracking, and agricultural field temperature-sensing networks [3]. Still another application 
consists of sensors that monitor conditions in the London underground tunnels and water 
systems [4]. 

Military applications may include remote sensing of nuclear, biological or 
chemical (NBC) agents, motion, video, etc. It is assumed that the sensor nodes are 
deployed in hostile territory and therefore delivery of the sensors may be via an aerial 
vehicle. Delivery in this manner would generally lead to a random spatial distribution of 
the sensors, yet the sensor nodes are close enough to form a local wireless sensor 
network. Once established, the issue becomes how to get the information from the 
wireless sensor network to the end user. The nodes of this network may act collectively 
as a group for such actions as distributed processing or forming a distributed antenna 
array. Vincent [5] and Chan [6] explored the use of beam forming to communicate 
directly with an overhead UAV as a method to relay information to the end user. 

In this thesis, we assume the deployed WSNs must communicate their 
information to the end user without access to a UAV. In the solution proposed here, 
where the WSNs are deployed to multiple locations nearby to one another, they form a 
terrestrial Over The Horizon (OTH) network by linking the individual WSNs. In order to 
form this link, the WSNs must also solve the additional task of locating one another. 
Each WSN acts collectively to form a smart antenna. A key assumption to this problem 
is that the distance between WSNs is greater than the individual range of a single sensor 
node. Thus the sensor nodes in the WSN must act cooperatively in order to have 
sufficient transmitting gain to reach the next WSN. 

This problem combines and draws from several different disciplines in order to 
propose a solution. Each of these topics will be discussed further in this thesis. In order 

to form the initial local network, the concepts of localization of the network are required. 
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Once the geometry is known, the sensor nodes of the WSN aet eolleetively to form a 
distributed smart antenna. A smart antenna generally refers to an array antenna that ean 
adapt its beam pattern using a eomputer proeessor [7]. Beeause the sensor nodes are 
seattered in a random manner, the eoneepts of random processes are used to apply them 
in the smart antenna composition. The concepts of spread spectrum and signal 
processing come into play in the determination of the direction of arrival for the beams 
eonneeting two WSNs. 


A. DESCRIPTION OF THE PROBLEM 

We begin with a deseription of the wireless sensor network, depleted in Figure 2. 
Consider a set of wireless sensor nodes, dispersed within an arbitrary boundary. The 
sensor nodes form an ad hoc network by establishing eommunications with one another. 
The boundary is defined as the range limit of the eommunieating sensor nodes within the 
ad hoc network. 


Sensor 


Sensor Node 
Connection 


WSN Cluster 
Boundary 



Figure 2. Wireless Sensor Network 


The individual sensor nodes may be randomly positioned within a boundary 
defining the WSN. The goal of this thesis is to explore methods for forming 
eommunieation li nk s between adjaeent WSNs. By establishing eommunieation links 
between adjaeent WSNs, we ean form an extended network made up of what has now 
beeome nodes of a larger Over The Horizon (OTH) network. The general ease is shown 
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in Figure 3. Once the OTH network is connected, the individual WSNs become part of a 
larger network that can relay sensed data to an end user. 



Sensed 

Phenomena 





WSN D 
End User 



WSN A 



WSN C 


WSNB 


Figure 3. Generalized Over The Florizon Network of Wireless Sensor Networks 


Next we consider a tactical military example where we wish to link WSNs to 
form an OTH network. Consider an Unmanned Aerial Vehicle (UAV) or other aircraft 
that disperses sensor networks to various locations in hostile territory, as envisioned in 
Figure 4. The goal is to connect the WSN A to B, B to C, etc. such that we eventually 
have a link from A to E. 
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Figure 4. Tactical Application of Linked Wireless Sensor Networks 

In previous work related to this scenario, each of the WSNs communicated 
directly to the UAV when it was available overhead [5]. In this work we seek to 
communicate not via a UAV, but by establishing links between already placed WSNs. 

In more general applications, a WSN need not be on the ground but can be 
attached to a building, bridge, structure or even a moving vehicle. The goal here is for 
the WSNs to find each other via beam forming and then establish a communications li nk 
at the highest rate possible. 


B. THESIS OBJECTIVE 

The goal of this thesis is to develop and validate methods to establish OTH 
communication links using wireless sensor networks by enabling each to act as a smart 
antenna array. 
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To begin forming the OTH eommunieation link, we examine methods to form a 
seareh beam and transmit this beam around the horizon in order to illuminate an adjaeent 
WSN. We will assume that eaeh WSN is eomposed of a uniformly distributed random 
array of sensor nodes. Using theory from periodie and aperiodie antenna arrays, we will 
develop methods to model and eharaeterize the random array in terms of gain and half 
power beamwidth as a funetion of the number of elements and the size of the random 
array. We will demonstrate using MATLAB a random planar array forming a beam and 
seanning the horizon. 

For the WSN being illuminated, the direetion of arrival (DOA) and loeation of the 
searehing array must be determined in order to form an optimum link. We seek to 
develop and simulate methods to determine the DOA and loeation information using a 
direet sequence spread spectrum based space division multiple access technique applied 
to the WSNs. 

In order for the WSN to act as a smart antenna, the local geometry of each WSN 
must be determined in order to define the sensor node positions and boundaries. Because 
this topic has been extensively reported in the literature, only an overview of applicable 
localization techniques is presented in Appendix A. To establish an optimum 
communications link, we will examine a pair of WSNs using a link budget analysis. We 
will determine the gain and number of sensor nodes required for this communication li nk . 
Finally, the overall OTH network is realized by establishing a series of communication 
links as described above. This area is present conceptually only, as the focus of this thesis 
is to investigate the establishment of the link between two adjacent WSNs. 


C. RELATED WORK 

The concept of beam forming using an array of randomly distributed sensor nodes 
has been explored previously in several works. Tong [8] studied random arrays in the 
context of randomly distributed micro-strip patch antennas placed on the sides of ships 
acting as a digital phased array radar aperture. Batson [9] studied wireless sensor 
networks as a means to collect signals intelligence by building adaptive arrays. Vincent 
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[5] considered a similar problem to the one we explore in this thesis, namely that of a 
WSN deployed in a remote location that uses adaptive beamforming to eommunicate to 
an overhead UAV. Vincent also explored energy burden sharing across the sensor 
network, a coneept that is coneeptually useful to an over the horizon network. Chan [6] 
studied adaptive beamforming for direetion of arrival ealeulation as well as the effeets of 
position error and element failure on the beam of a random array of elements. Godara 
[10] presents an extensive eompilation of teehniques for beamforming and direetion of 
arrival estimation, as well as a eomparison of the relative performanee of eaeh. While 
Hong, et al [11] use cooperative methods to establish a relay network of sensors, they use 
spatial diversity teehniques rather than beamforming. 

This thesis builds upon the concept of the random planar array of sensor nodes, 
but here the objeetive is to build an over the horizon network of WSNs. We also propose 
new ways for searehing, traeking and establishing a link between adjaeent WSNs. The 
communieation links between adjacent WSNs are established by using beam searching 
and direction of arrival methods based upon spread speetrum teehniques. 


D, THESIS ORGANIZATION 

This thesis is organized as follows. Chapter II provides a link budget analysis to 
determine the gain requirements for the transmitter and reeeiver given data rates, range 
and power considerations. Chapter III describes the use of beam-forming as a method to 
seareh for adjaeent WSNs. Simulation results for random planar arrays and seareh 
patterns are ineluded. Chapter IV deseribes the link establishment proeess and the 
determination of the direetion of arrival (DOA) using spread speetrum based spaee 
division multiple aeeess (SDMA) methods. Simulation results of a random planar array 
deteeting the arrival of a signal of interest as well as extraeting a sample data stream are 
presented. Chapter V presents the conclusions of this work and topics for further 
researeh. Finally, Appendix A summarizes the sensor node loealization methods and 
Appendix B includes MATLAB code used in this thesis. 
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II. COMMUNICATIONS LINK: LINK BUDGET 
CONSIDERATION 


The over the horizon (OTH) eommunieations system is essentially a network of 
networks linking the individual WSNs. Several issues must be resolved in order to link 
these WSNs of an OTH network. This thesis examines the link between just two of the 
WSNs, but for eompleteness, we discuss the bigger issues of the OTH network as a 
whole. Many of these issues are not too different from the issues discussed in Appendix 
A on localization. In this chapter, we begin with a conceptual look at the OTH network, 
then consider just two WSNs of the overall chain of WSNs that eventually form the OTH 
network. Within this view, we consider them as two points of a communications link and 
apply a link budget analysis. 


A, BUILDING THE OVER THE HORIZON NETWORK 

Each WSN of the OTH network must find the adjacent WSN to connect to. Once 
the initial link is established, say from A to B as in Figure 3, WSN B must now continue 
the link and connect to a third WSN, for example WSN C of Figure 3. While 
establishing the second link the WSN node must still maintain contact with the first 
WSN. A particular WSN may subdivide as a method to establish contact with more than 
one other WSN. The problem of subdividing an individual WSN may be addressed in 
several ways. Perhaps the WSN could physically partition the sensor nodes to address 
and connect to other nodes, e.g. the Eastern most nodes form an array to communicate to 
another WSN to the East, and the Western most sensor nodes form an array to 
communicate to the West, as depicted in Figure 5. In Figure 5, we show WSN B 
subdivided laterally, such that the red nodes form an array and communicate with WSN 
A and the green nodes form an array and communicate with WSN C. 
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Figure 5. WSN B Sensor Nodes Partitioned to Simultaneously Communieate 

With WSN A and WSN B 


Alternately, the WSN eould employ some subset of the sensors using a 
multiplexing teehnique based on time, code or space. Additionally, the frequency band 
could be split into different frequency channels for communicating to individual 
surrounding WSNs. This thesis does not address the details of subdividing the WSNs, but 
rather approaches this part of the overall process conceptually. 

In the following chapters, we will model the transmitting and receiving WSNs as 
smart antennas, taking advantage of the element localization information, the random 
nature of the element spatial distribution to form a transmit beam and a spread spectrum 
SDMA technique to determine the relative geometries between WSNs. 

Essentially, the WSN will now be functioning as individual antennas, and thus the 
communication link can broadly be described by Figure 6. This figure depicts a 
terrestrial communications link formed over some distance R between antennas A and B. 
Because of the nature of the WSNs, we can change some of the characteristics of the 
communication link, such as its effective aperture and the number of elements used to 
form the antenna. This chapter seeks to examine those parameters given the kn own 
parameters such as range, required data rate, path loss and available transmitter power. 
We begin with a discussion of the link budget equation and then perform an analysis 
using a commercial sensor mote and an example random planar array. 
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Figure 6. Simplified Terrestrial Radio Communications Link 


B. LINK BUDGET EQUATION 

The link budget allows us to perform a tradeoff analysis when we are given 
control over components of the link. In this section, we examine the individual terms of 
the link budget equation [12]; 



( 1 ) 


V o 


where Mds is the link margin (dB), Pt is the transmitter power (dBW), Gt is the gain of the 
transmitting array (dBi), G, is the receiving array gain (dBi), t is the energy per bit 

to noise power spectral density ratio (dBW/Hz), Rb is the bit rate of the communications 
between WSNs (dB-bit/sec), k is Boltzman’s constant (1.38x10'^^ JIK), is the 
effective system temperature (K°) generally given as T^=T^+ ,where Ta is the antenna 

noise temperature and Tg is the receiver noise temperature. The path loss, Ls (dB) is the 
loss of a signal over the distance from point A to point B, including multipath, fading, 
two ray bounce, weather, etc. and Lo includes “other losses” such as transmitter 
inefficiencies, line loss, polarization mismatch, etc. [12]. The link margin is the 
difference between what we need to establish the link and what we have available from 
the combined components of the communications link. We assume that each transmitter 
can vary its power output in discrete steps, but that all elements transmit at equal levels. 
The gain of the transmitting and receiving arrays refer to the gain of the random planar 
array of judiciously selected elements. The gain of a random array of transmitters will be 
described in Chapter III on random arrays and will be one of the key parameters we can 
vary. Here, we will assume the “other losses” are negligible. 
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1, Link Analysis 

Here we eonsider two WSNs deployed into separate but nearby regions on a flat 
environment, not unlike the scenario depicted in Figure 4. For example, let us consider 
as sensors the Crossbow MICAz, a commercially available sensor mote, to apply to this 
scenario. MICAz operates at 2.4 GHz and generates RF transmitter power from 40/^W to 
ImW. The nominal outdoor range for a single receiver/transmitter pair is given as 100m 
at the receive sensitivity of -120 dBW [13]. The maximum data rate for this mote is 250 
kbps. Using devices like these one may create a WSN containing 50 sensors arrayed over 
a 5 m X 5 m (40A x 40A) region. For simplicity, we will assume that both WSNs are 
the same size and contain the same number of elements and that the separation distance 
between the two WSNs is 500 m (4000 A). 

We will assume that all 50 of the receive sensor nodes are awake and in “receive” 
mode, although generally this need not be the case. We now evaluate the terms in the 
link budget equation for this example. The digital modulation scheme for the MICAz is 
BPSK or OQPSK per the IEEE 803.15.4 standard. As a first example, we assume that 
each transmitter is operating at the maximum power level of 1 mW which equates to -30 
dBW. At minimum power, the transmitters produce 40 juW or -54 dBW which we 
consider second. We assume that all the transmitters operate at the same power level. Eor 
the path loss term of the link budget, the contribution is not as clear. Before continuing 
with the link equation, we briefly examine the path loss for two sensor nodes positioned 
on the ground. 


a. Path Loss Models 

Many models exist for land mobile radios and wireless EAN (WLAN) 
systems and produce a wide range of values when applied to sensor nodes. A model to 
estimate path loss for sensor nodes at ground level is not readily available as most 
address antenna heights of approximately 1 m or above. Models that do include antenna 
height are relatively sensitive to the height of the transmit and receive antennas. Eor the 
MICAz sensor mote, which uses a one-half wavelength dipole (-6.25 cm) and sits on a 
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battery pack and circuit board ('-'3 cm), this gives us an approximate antenna height of 10 
cm. We now examine the most common models for path loss. 

The free space path loss at a distance d and wavelength X is given by [12]; 

4 = 101ogio — 

\ A J 

The free space path loss is not applicable to the scenario under 
consideration, but it serves as a baseline for comparison with the other models that will 
be presented below. In this case, the antennas are very close to the ground and therefore 
the path loss is expected to be greater. Green and Obaidat developed a path loss model 
for wireless LANs operating at 2.4 GHz that takes antenna height into account and was 
experimentally studied under conditions similar to ours. The Green and Obaidat path loss 
equation is given by [14]; 

Klan = 7.6 + 40 logio J - 20 logio hX (dB) 

where d is the separation distance and ht and hr are the heights of the transmit and receive 
antennas, respectively. 

Okumura [15] and Hata [16] produced landmark papers on propagation 
loss in land-mobile radios, but these models are for signals below 1500 MHz and 
generally for antenna heights of 1 m or more. The Green and Obaidat models 
demonstrated very good agreement with their experimental results, for antenna heights of 
1-2 m, but produce higher losses when used for antenna heights of 5-10 cm. 

A propagation model based on a flat Earth assumption, where the Earth is 
represented as a dielectric is another approach. Einnartz [17] provides the path loss at a 
distance d for the flat Earth as; 

L„p = --^ 2sm -—hX 

[AndfV yid ' X\ 

where X is the wavelength, and h and hr are the heights of the transmit and receive 
antennas, respectively. Because the antennas are so close to the ground and we are using 
a 2.4 GHz example this model would only apply in a very limited number of cases where 
the environment was free of scattering objects. 
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The calculated path losses for the 2.4 GHz signal at a separation of 500 m 
with both antennas at 10 cm are listed as follows: 


Free Space 

94 dB 

Hata Model 

121 dB 

Green/Obaidat Model 

134 dB 


Based upon the operating environment, the path loss for this system ranges from 94 to 
134 dB. Each of the models, except Free Space, is sensitive to antenna height. The Green 
and Obaidat model seems reasonable for use in this work given that their model was 
developed specifically for WLAN applications at 2.4 GHz and had good agreement with 
experimental results. However, this model used much greater antenna heights suggesting 
that using this model for sensor nodes on the ground may not produce accurate results. 

b. Link Budget Gain Calculations 

We wish to know what level of gain is required for the transmit and 
receive antenna arrays forming the communications link shown in Figure 6. We can solve 
the link equation to determine, with margin M, what gain the two WSN arrays must 
produce. Since we are given the receive power sensitivity Pr, of the MICAz, assumed to 
incorporate the data rate and required Equation (1) simplifies to: 

M,, = P,+G, + G^-P^-L^ ( 2 ) 

where M^b represents the difference between the actual received power and the minimum 
useable power, given by the receive sensitivity. Also note that the total power 
transmitted by the array of Nt transmitting elements is N,P,. 

Using Equation (2) and combining the terms from the link analysis and the 
Green/Obaidat path loss estimate, we have for the 1 mW (-30dBW) transmitter: 

M = (G, + G^) + N,+P,-P^-L^ 

M = (G, + GJ + A, - 30 - (-120) - 134) 

(G,+G^) = M-N, + 44 

Thus the combined gains of the transmitting array and receiving arrays 
must meet a margin balanced by the total transmitted power less a factor for path loss. 
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Excluding the free space model, the range of path loss values is on the order of 15 dB, 
suggesting that we should ehoose this as our margin. 

When the transmitters are operating at minimum power (-54 dBW), from 
Equation (2) the required gain is: 

(G,+GJ = M-7V,+68 

In this ease, it is clear that we must now either decrease our acceptable margin or increase 
the number of transmitters or both in order to form the communications link. Because 
transmission consumes much more energy than reception, it is envisioned that we will 
use many more reeeive nodes than transmit nodes to meet the gain requirement. Eor 
example we may use 15 transmit and 30 reeeive nodes. In Chapter III, we will be able to 
determine these numbers. 

Clearly the path loss dominates the overall margin requirement. One 
method to reduee the path loss is to elevate the sensor nodes. This increases the antenna 
height and rapidly reduees the path loss, thus leading to greater range. Typically, 
doubling the antenna height inereases the reeeived power by a factor of 4 [17]. Eor 
example, the sensor network could be plaeed on the side of a building or rooftop, 
elevated on top of a vehiele, attached to lamp posts, etc. 

In summary, we applied the link budget equation to a link formed between 
two WSNs and determined the required combined gain from the reeeive and transmit 
antennas. In the next ehapter we examine how to ealeulate this required gain as well as 
some other parameters of the antenna we will form from the WSN. 
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III. THE SEARCH FOR OTHER WIRELESS SENSOR 

NETWORKS 


In this chapter, we approaeh the communications link from the viewpoint of a 
searching WSN in the eventual OTH network. For example, this could be WSN A of 
Figure 3, which has sensed a phenomenon and seeks to pass this information along to the 
end user. This WSN has the task of transmitting a search beam in order to find adjacent 
WSNs. Like the rotating beam of a lighthouse, the WSN will steer its beam in a rotating 
manner around the horizon by using the loeal network of sensor nodes as an aperiodic 
antenna array. 

In order to address the method, we first discuss the concepts of periodic arrays. 
Then we extend these concepts to the random array, a type of aperiodic array. Next we 
explore the charaeteristies of the two dimensional random array, using MATLAB to 
calculate gain, beamwidth and sidelobe levels. Finally, we discuss the issues of the 
seareh between WSNs, using random planar arrays. 


A, BEAM FORMING WITH PLANAR ARRAYS OF RANDOMLY PLACED 

SENSOR NODES 

To describe how we will form a beam using a group of randomly placed sensor 
nodes, we begin by reviewing the theory of the uniform linear array (ULA) and the 
uniform planar array. From this, we make the extension to planar arrays containing 
randomly placed elements. In each case, we assume that the sensor antenna elements are 
half wave dipoles oriented vertically to the plane. 

1. Geometric Arrays 

The far field pattern for an array of geometrically arranged antennas is caleulated 
using pattern multiplication; that is, the far field radiation pattern F{6) can be expressed 
as the product of the element pattern £^{6) and the array factor AF(0) [ 18 ]: 

F(e)=E„{e)xAF{e) 
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For a z -directed, half wavelength dipole, the normalized element factor is [18]; 


cos 


EM-. 


K ^ 

COS0 
2 ) 


sin0 


( 3 ) 


and is assumed to be the element in all the gain pattern calculations that follow. The 
array factor characterizes the radiation pattern based upon the geometry of the antenna 
elements. We now discuss the one-dimensional uniform linear array. 


a. Uniform Linear Array (ULA) 

Figure 7 shows the geometry of a uniform linear array and its coordinate 
system, and observation direction R(0). The angle 6 is measured from the z axis. 


Z 



Figure 7. Uniform Linear Array Geometry and Coordinate System 

For an equally spaced set of N isotropic radiators with wave number 
k — lKlX, spaced by a distance d along the x axis, the array factor may be expressed as: 
[18] 

N-\ 

= (4) 

«=0 
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where is the amplitude of the drive of the eurrent element, and is the phase 

eontribution of element. Note that the array factor is rotationally symmetric around 
the x-axis. Figure 8 illustrates how the observation angle 6 and the spacing d affect the 
phase term of the array factor. For example, when 6 is zero, all the elements have the 
same phase because sin(0) is zero; when sin(0) = 1, the phase contribution from each 
element is equal to kd. From this, we see that the spacing, observation angle and wave 
number determine the phase of the far field array factor. 



Element 


Figure 8. Phase Difference is determined by 0 and spacing d 


We can simplify Equation (4) by letting 

i// = kdsmd 

and letting all the current amplitudes be equal, i.e., 

A) = A = - = ^v-i 

Thus, Equation (4) becomes; 




n=0 


( 5 ) 


For elements, and the beam pointing in the 6=0 direction, the array factor reduces to 

AF{y/) = A,N 

Thus, the maximum array factor is proportional to N. 
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It can be shown [18] that Equation (5) can be expressed as the more 
common expression for the normalized array factor /, or pattern factor of a uniformly 


excited, equally spaced linear array, centered at the origin: 

smjNwn) 

^ ’ A^sin(v//2) 


( 6 ) 


A close examination of Equation (6) reveals some important characteristics of array 
antennas [18]: 


• As the length of the array {Nkd) increases, the width of the main lobe 
decreases, i.e., the half power beamwidth (HPBW) becomes smaller. 


• As A increases, the side lobe levels decrease. The Side Eobe Eevel (SEE) 
is expressed as the ratio of the maximum value of the largest side lobe to 
the maximum value of the main lobe. This is usually described in the 
literature as “dB below the main lobe.” 


A characteristic of periodic arrays is revealed when the spacing between 
elements increases beyond a certain level, and energy from the array is directed into 
unintended directions. These unintended beams from the array essentially reduce the gain 
of the array along the desired direction by putting energy into unwanted directions. 
These unintended beams are called grating lobes and are the subject of the next section. 


b. Grating Lobes 

Grating lobes occur when element spacing is greater than one-half 
wavelength [18]. These lobes are often of equal intensity to the main lobe(s) of the array 
and are due to the periodicity of the element spacing. Grating lobes do not occur for 
arrays with spacing less than one-half wavelength; thus, nearly all periodic arrays 
conform to this spacing constraint [18]. In order to demonstrate the occurrence of grating 
lobes, Eigure 9 shows the pattern factor for an eight element UEA, with a spacing of 0.5, 
1.0 and 1.25 wavelengths between elements. Only the first quadrant is shown, as the 
pattern is symmetric about 0°. The blue line shows the pattern factor for an array with 
one-half wavelength spacing and no grating lobes. The red and green lines show the 
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pattern factor for arrays with a spacing of 1.0 and 1.25 wavelengths and the resultant 
grating lobes. Notice that for the one-half wavelength spaced array, shown in blue, most 
of the energy of the array is concentrated in the main beam, while for the larger spacing 
the energy is dispersed into the grating lobes. 



Figure 9. Pattern Factor for an 8 Element ULA Showing Grating Lobes as a Result of 

Spacing Greater Than One-Half Wavelength 

Thus the designer is presented with a dilemma - in order to have a small 
beamwidth, the ULA must be physically large, but increasing the length will necessarily 
result in either a higher complexity due to the greater number of elements or, if spread 
out too much to keep the number of elements down, grating lobes. It will be shown later 
that since random arrays do not have a periodicity, grating lobes generally do not occur. 

We now extend the concepts of the uniform linear array to form a two- 
dimensional periodic array. 

c. Two-Dimensional Periodic Arrays 

Ligure 10 shows the geometry for a two-dimensional periodic array. 
Elements are geometrically arranged in the x-y plane, and the far field observation 
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direction is shown by R{6,(j)). For this array, the principle of pattern multiplication is 
also applied. Here we can write the far field pattern F(0,(j)) = F^(0,(j)) x Fy{9 ,^), where 
is the pattern of a linear array along the x axis, and F^{0,(p) is the pattern of a 
linear array along they axis [19]. 


Z 



Figure 10. Uniform Two Dimensional Array Geometry and Coordinate System 


The array factor for a two dimensional periodic array with M elements of 
spacing dx in the x direction and N elements of spacing dy in the y direction is given by 
[ 20 ]: 


AF(0,0) = ££a„ 


^ jjt (m-1)(rf, sin 0 cos i/)+ft )^ (n -1 sin e sin i/i+/J,, ) 


m=l n = \ 


( 7 ) 


Combining the exponential terms yields: 

M N 

. jk[{m-l)(d^ sin0cos0+/3^. )+{n-l){dy sindsin^+^y )] 

m=l «=1 


AF(0,0) = XX 


( 8 ) 


where Amn is the current distribution at the element. The array can be steered to an 

arbitrary direction (6^,(1)^) via setting [20]: 
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= -kd sin0^ cos0^ and = -kd sin0^ sin0^ (9) 

When combined with a specific radiation element, the total radiation 
pattern is found by multiplying the element pattern £'q( 0,^) with the two dimensional 
array factor AF(d,(j)) [18]; 

F{e,(^) = E,{G,(l))xAF{e,(p) (10) 

The element pattern is the radiation pattern of a single element in 6 and <p. In our case, the 
radiation pattern for a half wavelength dipole given is rotationally symmetric in the angle 
(f); therefore. Equation (3) can be used in Equation (10) as {9,(p). 

It can be seen that with half-wavelength spacing, even a modest 7x7 
wavelength array would contain nearly 50 elements. It is possible to thin out these arrays 
yet still have control of the beamwidth, gain and grating lobes. One can thin the arrays 
by placing the elements randomly within a given area at a much lower density than that 
of a periodic array. The next section develops these concepts. 

2. Random Arrays 

A random array is an array antenna whose elements are no longer defined by a 
fixed geometric spacing but rather the element locations are now random variables. 
Eigure 11 depicts the periodic linear array, the random linear array and the random planar 
array. A one-dimensional periodic linear array has a fixed inter-element spacing, while 
the element spacing in a random linear array is a random variable. The elements are 
randomly placed along both the x and y-axes in a random planar array, i.e., the inter¬ 
element spacing is random along both axes. 
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Figure 11. Illustration of Periodic Linear Array, Random Linear Array and Random 

Planar Array Showing Example Element Distribution 

The array factor for a random linear array is similar to the array factor for a 
periodic array, repeated here for convenience; 

N-l 

«-0 

except that now the element spacing term in the exponential is a random variable. Eor 
the purposes of performance comparison, the periodic array is termed the design array 
and its array factor is the design array factor AFq{0) . The phase contribution from each 

element in a random linear array is no longer a deterministic term, i.e., nd ^, but rather 
from the actual x coordinate of the n^ element, x,„ and thus the phase is also a random 
variable. The array factor for linear random array becomes: 

iV-l 

= (11) 

rt—0 

Because the element location is now a random variable the array factor is also a 
random variable. It is customary to represent the array factor of a random array in terms 
of the ensemble average, referred to as the average array factor, given by: 

AF{e) = E[^AF{e)~] 

where £’[•] is the expectation operator. 
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As two dimensional, or planar arrays, are eentral to the foeus or this work, we 
extend further the random linear array to the random planar array. In the same manner as 
with the linear arrays, the inter-element spacing is no longer deterministic, but random in 
both coordinates. In the previous discussion about periodic planar arrays, the array size 
contained N x M total elements, but now we express the total number of elements as 
simply N. The double sum in Equation (8) is merely a convenient way to express the 
sum over the two dimensions. Here, it is a single sum because we explicitly step through 
the coordinates of each element in a random planar array. Thus, for N total elements, the 
array factor becomes: 


AF(0,0)=XA„e 


sin0cos0+^^)+(_y,, sin6sin0+/3 )] 


n=0 


( 12 ) 


Substituting 


= -kx„ sin0„ cos0„ and = -ky^ sin0^ sin0^ 


(13) 


into Equation (12) gives the array factor for the N element random planar array: 


AF(0,0)=XV 


Jk[{x^ sin0cos0-x„ sin 6^ cos0^j)+(}’« sin 0sin0-_y„ sinQ^ sin0y)] 


(14) 


n=0 


where {6^,^^) determine the beam pointing angles. This equation is used to determine 
the array factor in all subsequent calculations involving random planar arrays. 


Random arrays exhibit some special characteristics that differentiate them 
from periodic arrays. Some details of random arrays are presented in the following 
subsections. We will focus our attention on the array performance characteristics of 
grating lobes, the average array factor, main lobe and beamwidth, average side lobe level 
and peak side lobe level. 


a. Grating Lobes 

Grating lobes are a result of the periodic spacing of the elements of the 
array; there is a coherent buildup of elemental signals at angles other than the steering 
angle. It should be no surprise then that grating lobes, an artifact of the periodicity of 
uniform arrays, are significantly reduced in random arrays because of the random nature 
of the elements [19]. 
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b. Main Lobe and Beamwidth 

The shape of the main lobe (maximum gain and HPBW) and the shape of 
nearby side lobes are not signifieantly affeeted by the random nature of the elements [19]. 
This can be shown by examination of the expression for the average power pattern of a 
linear random array: 

( 15 ) 

V iS J rs 

where |/o(0)f is the design array power pattern and N is the number of elements. More 
precisely, near the main beam where the normalized power pattern of the design array is 
approximately one ^l), the main lobe character of the random array (maximum 

gain, half power beam width) is approximately that of the periodic array. Furthermore, 
the variance of the array factor is given by [19]: 

where /o(0)is the design array factor and N is the number of elements. As the number of 
elements N grows, in the limit the variance approaches zero: 

lim (jI =0 

This simply states that the larger the number of elements, the smaller the variance, and 
thus the closer the power pattern of the random array is to the design array power pattern. 

The broadside half power beam width (HPBW) of a linear geometric array 
is given as [18]: 

X X 

HPBW = 0M6-rad = 50.- (16) 

L L 

We wish to compare how the HPBW of a random array compares to the HPBW of a 
design array of the same length. 

Figure 12 graphically shows how Equations (15) and (16) compare to the 
power pattern and HPBW of a linear array - the design array; shown here for arrays of 5, 
8 and 20 elements with fixed inter element spacing. As expected, the power pattern of 
the random array is almost the same as the design array in the vicinity of the main beam. 
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but forms a “pedestal” of magnitude 1/N as evidenced in Equation (15). The random 
nature of the array will cause the average side lobe level to be equal to the pedestal. We 
also see that in the vicinity of the main lobe, the HPBW is approximately equal to the 
design array. The approximation improves as the number of elements increases. 




- N= 8, (1(*= 0.5 

-Random Array Ave PP 


0 (deg) 


-20 - 



Figure 12. Comparison of the Broadside Power Pattern of 5, 8 and 20 Element Einear 
Array and Average Broadside Power Pattern of the Same Sized 
Random Einear Array 


The results of Figure 12 are important because they give us an envelope 
for the expected power level of the random arrays to be developed to build the proposed 
communication link between WSNs. We also see that as the number of elements 
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increases, the HPBW of the random array approaehes that of the linear array. This is 
important as we ean then easily estimate the HPBW of a random array based on using a 
simple formula for an N element linear array. Therefore, the half power beam width of 
the main lobe, for either a random array or a uniform linear array with equal spaeing and 
eurrent is given by Equation (16). 

c. Average Side Lobe Level 

While the main lobe oeeurs in the direetion the array has been steered, the 
side lobes are those that oecur in all other direetions. The arithmetie mean of the peaks of 
these side lobes is the average side lobe level. In eontrast to the structure of the main 
lobe, the distant side lobe strueture is affeeted by randomness of the element spaeing. 
This can be seen from the array faetor in Equation (14), where Xn and ya are random 
variables, thus the phase is a random variable. Off of the main beam, where the power 

pattern of the design array is mueh less than one (|/o(0)|^ «1)? we see that from Equation 
(15) the power pattern of the random array is dominated by the J/N term. Sinee the 
maximum of the power pattern of a random array is equal to that of the linear array, and 
the average side lobe level of a random array is 1/N, we ean say that the power in the 
average side lobe level is equal to iW below the main lobe [19]. In Eigure 12 we can see 
this graphieally. All plots are normalized to 0 dB. In Eigure 12 we see the average 
sidelobe level for the 20 element array is approximately 1/N = 101ogjo(l / 20) = -13 dB 
down from the main lobe. 

d. Peak Side Lobe Level 

While the maximum of all side lobes of a periodie array is generally 
termed the side lobe level [18], for the ease of the random array all of the sidelobes are 
random in nature and average to 1/A. The peak side lobe level refers to the maximum 
peak sidelobe value relative to the average side lobe level [19]. The topie is not easily 
summarized; however, a rule of thumb for the peak side lobe level is that it is unlikely to 
exceed the average side lobe level by 10 dB [19]. 
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In conventional periodic arrays, an increase in performance can be 
achieved by using a non-uniform current distribution across the elements of an array. For 
a given number of elements, by tapering the current amplitude near the edges of an array, 
the side lobe level of an array can be reduced and the main beam width increased [18]. 
For random arrays, however, this is not the case; a tapered current distribution has only a 
slight impact on the side lobe level [19]. Therefore, in our simulations we will use 
uniform current amplitude at all nodes of the antenna array. 

Having discussed the characteristics of periodic and random arrays, we now 
proceed to demonstrate this theory using a MATLAB simulation of the random planar 
array. 


B, MATLAB SIMULATION RESULTS - RANDOM PLANAR ARRAY 

In this section, the periodic planar array with fixed element spacing is extended to 
uniformly distributed, randomly placed elements over a planar area. The elements are z 
directed half wave dipoles. The array size is characterized in terms of wavelengths and 
the number of elements N, which represents the total number of elements. This is in 
contrast to periodic planar arrays which commonly list the number of elements along 
each axis, e.g., Nx= 5 and Ny= 5 for a total of 25 elements. Here, we simply use N= 25 
to indicate that there are 25 elements within the array. 

Elements in the random planar arrays are restricted to being no closer than 1/10^*' 
of a wavelength. This is done for two reasons, to better model a physical sensor node as 
it is unlikely that one would sit on top of another and also to reduce the influence of 
mutual coupling between antenna elements. Although at a separation distance of 1/10^*' 
of a wavelength, there is still a contribution from mutual coupling [18]; in this work it is 
ignored. 

In Figure 14, we show the 3 dimensional gain pattern of a random planar array of 
25 vertically oriented dipoles with the elements constrained within a 3Ax 3A area. This 
plot is presented to demonstrate the pattern of the random planar array in the upper 
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hemisphere, with the array beam steered to 9= 90° and (p = 0°, here oriented to the lower 
left. Higher gain is shown in bright red, while lower gain is shown in purple and blue. In 
this thesis, we are interested in forming the terrestrially bound eommunication link and 
thus we will foeus on the portion of the pattern that lies in the x-y plane. 



Figure 13. Three Dimensional Gain Pattern of 25 Element 3Ax 3A 
Random Planar Array Steered to 6= 90° and 0=0°. 


In the following seetions, through simulation we examine the relationships 
between the physieal array size and the resulting half power beamwidth; the gain of a 
fixed size array and the number of elements eontained within it; and the gain of the array 
with a fixed number of elements versus the array size. These results will be useful for 
developing the proposed eommunieations link between WSNs in later chapters. Note 
that since the WSNs are all assumed to be on the ground, the search will be along the 
horizon, or in our simulation’s coordinate system (see Figure 10) we have 6 = 90°, and <p 
scanning the x-y plane. For the actual sensors (see Figure 1) the transmit and receive 
antennas are vertically oriented, thus we only simulate the vertical polarization. 
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Additionally, we assume here that the radiation effieiency, r]is one so that the terms 
directivity and gain are used interchangeably. 

1. Planar Array Size versus Half Power Beam Width 

Given a random planar array of half-wave dipoles parallel to the z-axis, we wish 
to find the half power beam width in the direction the array has been steered. The 
broadside half power beam width of an array with wavelength A and length L is: 

A 

HPBW = 50.8° —. This also applies to a random linear array of length L [19]. We now 
L 

consider the uniform random planar array, steered to 6o = 90° and <po = 0°. For 
observation angles in the x-y plane, 6= 90°. Here we are interested in observations angles 
near the main beam, i.e. where <p is small and coscj) = 1. Thus Equation (14) reduces to: 

N-l 

AF(e = 90°,0) = ^ 

n-0 

where x„ and range from 0 to T and all are equivalent. When 0 = 0° this is then the 
broadside of the array where the gain is maximum. For observations where <p is small, the 
length L in the y dimension determines the HPBW. Because cos0 changes much more 
slowly than sin0, the contribution due to the x component is relatively small. Thus the x 
component of the element has a very minor influence on the HPBW in the broadside 
observation. Therefore the expression for the HPBW of a linear random array is 
approximately valid for this case of the random planar array. Because the distribution of 
elements within the Lx L area is random, this makes the choice of coordinate system 
arbitrary. If the pattern is constrained to the x-y plane, the expression for HPBW is valid 
for an arbitrary steering angle in the x-y plane. For larger values of 0, the array factor is 
dominated by the random character of the elements as demonstrated in Figure (12). 
Having shown that the HPBW of an N element random array can be estimated using 
Equation (16), we now validate this through MATFAB simulation of the random planar 
array. In this simulation, we model the planar array by arranging elements with a 
uniformly distributed random spacing. In this simulation, we held the number of 
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elements constant at 20. We then varied the length of the sides of the array, in terms of 
wavelength and calculated the HPBW at (9 = 90° and <p= 0° for the vertical polarization. 
The coordinate system is the same as in Figure 10. For each array size, we generated 50 
realizations of the array of random elements and then averaged the HPBW values from 
each to form a single data point. At each array size. Equation (16) was evaluated and 
plotted as the theory. The results of the simulation are presented in Figure 14. 
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Figure 14. Broadside Half Power Beam Width for 20 Element Random Planar Array and 
Theoretical Half Power Beam Width Calculated for Increasing Array Size 


We can see that as the array size increases, the HPBW becomes smaller. The 
results show good agreement with the theory, indeed the analytical expression seems to 
form an upper bound. Prom this result, we can use the array size to vary the beamwidth 
when the WSN is either searching for other WSNs where we want a relatively wide 
beam, or communicating directly with another WSN where we want a high gain narrow 


32 













beam. It should be noted that as the array size grows the spaeing between elements will 
become large relative to a wavelength, especially when the array size is 20A x 20A . 

2. Planar Array Maximum Gain versus Number of Elements and 
Average Side Lobe Level 

In this section we examine how the maximum gain of a planar array is affected by 
increasing the number of elements for a fixed size array. In Section III.A.2, we observed 
that the gain of an N element random linear array approximates that of an N element 
linear (design) array; additionally, the average side lobe level is 1/A below the maximum 
gain and the peak side lobe level is no more than 10 dB above the average side lobe level. 
This simulation will validate these observations for the random planar array. 

In Figure 15, we examine the maximum gain of the N element random planar 
array of half wave dipoles versus the number of elements. The physical size of the array 
is held constant at 5A x 5A . The maximum gain is calculated for the steering angle 6o = 
90° and 0o = 0°, and is plotted in blue. Also plotted are the average side lobe level in red 
and peak side lobe level in cyan. The average side lobe level was computed by averaging 
all the peaks of the side lobes for each given number of elements. For each number of 
elements, 75 realizations of the random array were computed in order to form an average 
value and the datum plotted. The green line representing the theoretical average SLL is 
from Section III.A.2.c and is estimated as I/A below the maximum gain. The purple line 
is the theoretical gain for a linear array of A dipoles and is lOlogio(A) plus the gain for a 
half wavelength dipole of 2 dB [18]. The peak side lobe level in cyan is presented for 
reference. As discussed in Section II.A.2.d above, the peak SLL is less than 10 dB above 
the average SLL. The cyan line shows this and is clearly in agreement with theory. 
More exact, but significantly more complicated predictions for the peak SLL are 
presented in Chapter 8 of [19]. 


33 




Figure 15. Maximum Gain, Average Side Lobe Level, and Peak Side Lobe Level as a 
function of the Number of Elements in a Random Planar Array of Size 5X x 5X. 

Theoretical Values Are Shown for the Maximum Gain and Average SLL. 

In Figure 15 we see that the maximum gain is increasing with an increasing 
number of elements in the array and that as the number of elements increases, the average 
SLL is decreasing farther and farther below the maximum gain. The peak side lobe level 
is also decreasing as the number of elements increases. This is in agreement with theory. 

Note that the maximum gain comes from the main lobe of the array, and is 
therefore unaffected by the random nature of the array, but rather is proportional to N 
[18]. This is in agreement with the plot. The difference in gain for 30 or more elements 
is less than or equal to 1 dB and the simulation and theoretical results converge as N 
increases. Overall, the theoretical and simulation curves are clearly in agreement. 
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From the theory of the random array, the average side lobe level should be less 
than the main lobe by a faetor of UN. This is depieted by the red and green lines in Figure 
15. The theoretieal and simulation results again show exeellent agreement. For example, 
at = 20, the average SLL should be 13 dB below the peak. Here the peak is 
approximately 17.0 dB while the average is approximately 4.7, yielding a differenee of 
12.3 dB. At A^= 50, the differenee is 16.3 dB versus a ealeulated 17 dB. 

These results show that we ean vary the amount of gain for a fixed size array by 
varying the number of elements we choose to include. Thus, for a particular gain level 
required to form the proposed communications link, we may save battery power by only 
using the minimum required number of sensor nodes to transmit the message while 
simultaneously extending the overall lifetime of the WSN. 

3, Planar Array: Gain versus Array Size for a Fixed Number of 
Elements 

In a periodic array, gain increases with array size because elements are uniformly 
spaced. In a WSN, we may wish to achieve a particular gain level without increasing the 
number of elements. In this simulation, we demonstrate that we can thin a random planar 
array and maintain a constant gain. We fix the number of elements in the array at 20, 
then increase the size of the array from a lA x lA to a 20A x 20A . At each array size, we 
compute the maximum gain, average SLL and peak SLL as was done in Figure 15. 

Figure 16 shows the results of this simulation. We note that for a periodic array, 
the gain is actually due to an increase in the number of elements. Here we see that for a 
fixed number of elements (20), the gain shown in blue is relatively constant over a 
varying physical size. Since the number of elements is held constant and the physical size 
is increasing, we can see the effect of reducing the element density, or alternately, 
thinning the array. Also note that the average side lobe level shown in red is 
approximately 13 dB below the main lobe gain, which is given by Equation (15) and 
shown in green. As expected, the peak SLL does not exceed the average SLL by more 
than 10 dB. This is shown in cyan. 
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Figure 16. Maximum Gain, Average and Peak Side Lobe Levels as a Function of 
Random Planar Array Size for Fixed Number of Elements = 20. Theoretical 
Values Are Shown for the Maximum Gain and Average SLL. 

As in Figure 15, the measured maximum gain is in agreement with theory; 
however, here it stabilizes approximately 2 dB above the theoretical level for a 20 
element array of dipoles of 15 dB, shown in purple. The simulation appears to break 
down for array sizes less than 4 wavelengths, yet it quickly converges to a constant gain. 

The independence of gain versus array size is an important result. For the WSN, 
we may require a particular gain level as determined by the number of elements N 
composing the array. From the results of Figure 16, we may choose any N elements as 
long as their distribution is still random within the bounded area [19]. 

4, Summary of Results 

Figures 14, 15 and 16 demonstrate some important results that we can employ 
when building the proposed communications link between WSNs. The particular WSNs 
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are assumed to be fixed in loeation, but the inter-element spaeing follows a uniform 
random distribution within a defined boundary. In order to link the WSNs, we wish to 
eontrol the width, gain and direetion of the array beam. We have shown in Equation (14) 
that we ean steer the beam in an arbitrary direetion. In Figures 14, 15 and 16 we have 
demonstrated eontrol over the beamwidth and gain. 

To achieve a given beamwidth, we choose a physical size, likely a sub-section of 
the WSN under consideration. The physical size controls the HPBW. Figure 14 and 
Equation (16) demonstrate this behavior. Within this physical area, we have control over 
the gain as demonstrated by Figure 15 by choosing the number of elements employed in 
the array. 

Alternately, for a given required gain level. Figure 16 shows that we can thin the 
array (or reduce density) dramatically without significant loss of gain or the effects of 
grating lobes. It does not matter which particular sensor nodes within the WSN are used 
as long as, on average, they are randomly distributed within the physical area chosen. 
This is important as we may wish to distribute the burden of transmitting, and thus 
consuming energy among the nodes within the WSN. 

To continue from Chapter II, in the first example using a pair of 50 element 
random planar arrays, with transmitters operating at full power (ImW each) and a margin 
of 15 dB, we determined the combined gain to be (G, -l- G^) = 15 - A, -l- 44 dB. From 
Figure 15, we see that by choosing A, =25, we have a transmitter array gain G? of 
approximately 17.5 dB. Therefore, we require a receiver array gain Gr of approximately 
16.5 dB. This can be achieved with 17 receive elements. In the second link analysis 
where the transmitters are operating a minimum power, the required combined gain must 
satisfy (G, -I- G^) = 15 - A, -I- 68 dB. Again, using Figure 15, we see that using 45 
transmitter elements gives us a transmitter array gain G? of 19 dB. The required receiver 
array gain Gr is approximately 19 dB. This can be achieved with approximately 45 
receive elements. 

Also consider a 20A x 20A area with 240 elements contained within the 
boundary. Suppose that the gain requirement is such that we only need 15 elements. We 
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are still free to ehoose a beamwidth based upon the eommunication link needs. Suppose 
that a beamwidth of 10° is required, which by Equation (16) implies a 5Ax5A sized 
array. Clearly, the WSN may subdivide itself into 5 A x 5A segments that contain at least 
15 elements. The sub-sections could then time share the burden of transmitting. Vincent 
in [5], provides methods for power conservation and burden sharing of the WSN. 

Armed with these results for determining the beamwidth and gain, we are now in 
a position to apply them to simulate a WSN searching for other WSNs. 


C. SEARCH METHOD 

In order for one WSN to find another, a method must be employed to 
systematically search the horizon. In this thesis we have chosen to employ a “lighthouse” 
approach in that we form a narrow beam and transmit in a given direction, then steer the 
beam across the horizon in search of an adjacent WSN. This section details some aspects 
of this method, such as beam steering, allowing for response time, power consumption, 
range capability of the array, calculation of the relative range and bearing between 
WSNs. The content of the initiating message is also discussed here. 

1. Beam Steering 

In this chapter, we have established that a random planar array can be used to 
form a beam, and this beam can be steered to a given direction. Equations (13) and (14) 
allow us to steer the array and the beamwidth is given by Equation (16). We know that 
beamwidth is related to the size of the chosen array (the longer the baseline, the thinner 
our beamwidth and vice versa), and that gain is related to the number of elements used. 
Eor example, to scan the horizon we may choose a half power beamwidth of 30°, transmit 
a signal in search of another WSN, then steer the beam ahead 30° and transmit again, etc., 
until we have moved around the entire horizon. 
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Look - Wait - Rotate Beam 


This search method is such that the searching WSN will pick a start direction, 
transmit a signal in that direction and wait for a response. If no response is received, then 
the searching WSN will advance the beam and “look” in the next direction. 

Since each node in the WSN acts as one element of the smart antenna, it has a 
particular phase or weight associated with it to steer the beam in each look direction. The 
weights can be stored in a table and sequenced through in order to advance the scan 
direction. These steering weights are the P^and of Equation (13). Since the 

geometry of the WSN is fixed, the steering weights associated with each sensor node (or 
transmit element) do not change and therefore do not need to be continuously 
recalculated. All nodes within the WSN transmit with the same amplitude. As noted in 
Section III.A.2.d on random array theory, tapering the current amplitude weights An of 
Equation (12) with respect to location in the array has little affect on the performance of 
the array. Since all elements transmit with the same amplitude, no complex current 
amplitude weighting scheme need be designed, simplifying the beam steering. 

The technique of rotating the beam and only looking in one direction at a time 
was chosen to avoid having several WSNs respond at once as might occur if a single 
omnidirectional beacon were used. If a single beacon signal were sent omnidirectionally, 
the searching WSN would need to determine the direction of arrival from each of the 
responding signals, possibly at the same time. 

3. Power Conservation 

Power conservation is not a major focus in this thesis; however, the results of this 
work are meant to complement other research on energy conservation in WSNs [5], [9] , 
[ 21 ]. 

When initially searching for another WSN, the array does not know the distance 
to the other array, and therefore must expect the longest distance possible. This implies 
using all elements available within a subset of a given area to form the initial array, to 
maximize gain. Since we initially want a relatively wide beam, for fewer total beams in 
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the search and thus fewer transmissions, these constraints limit the physical size of the 
employed sensor transmitters. Thus a tradeoff will be required to optimize power 
consumed to find another WSN and establish the link. 


4, Range and Bearing 

It is assumed that one goal of the WSN is to extend its lifetime as long as feasible. 
Since the sensor nodes are battery powered, and using the radio transmitter consumes a 
relatively large amount of power [5], we wish to minimize the number of transmitters 
used in forming the communications link. We wish to only use the minimum number of 
transmitters required for a given range. Since a higher number of nodes equates to a 
higher transmission gain Gt, from the Friis equation, we can determine the amount of 
gain required for a given distance J [12]; 

( ^ Y 

n= — P,G,G^ (17) 

\‘Xnd y 

where Pr is the received power, Gr is the receiver gain, X is the wavelength and d is the 
range. From Equation (17), we see that a higher transmitter gain produces a higher 
received power Pr, thus for a given receiver gain a pair of WSNs will be able to 
communicate at a longer range if more nodes are used. 


It is assumed that within the initial exchange of information between the 
searching WSN and receiving WSN are the respective latitude and longitudes. From this 
pair of lat/long coordinates, each WSN must calculate a bearing and distance to the other. 
For a given pair of latitude and longitude coordinates {laty,longy) and [lat 2 ,long 2 ), the 

distance between them is given by the Haversine formulas [22]. Here, the WSNs are 
assumed to be on the order of a kilometer apart and the fiat Earth approximation to the 
Haversine formulas is applicable [22]. The separation distance D (in km) and bearing 
O (in radians) are calculated using Equations (18) and (19): 



O = tan ‘ 


^ sm{Along)cos[lat2) 

^ cos{lat ^) sm(lat 2 ) - sin(/ah) cos{lat 2 ) cos [Along ) 


(18) 

(19) 
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where Along = long 2 - long^. A simple java-seript of the 

Haversine formulas for the flat Earth assumption is available at [23]. Once the range and 
bearing are known, the WSNs can calculate the required gain to establish the link. This 
gain calculation was the subject of Chapter II. For the sensor nodes, presumably 
deployed on the ground, the antenna height hant is approximately 0.1 m. The radio 
horizon for this node is J = 4.12^//z^ [24] or 1.3 km. Thus the flat Earth approximation 
used in Equations (18) and (19) for this short distance seems applicable. 

We now demonstrate a specific simulation of a WSN forming a beam and 
scanning it across the horizon. 


D, MATLAB SIMULATION RESULTS - SCANNING BEAM 

We now wish to combine the concepts of this chapter on random planar arrays 
with the idea of a searching WSN that transmits a beam in a rotating fashion across the 
horizon. We begin by choosing a half power beamwidth of approximately 20 degrees, 
from Equation (16), and set the array size at 2.5Ax2.5A. In this example we do not 
have a specific gain requirement to meet, so we pick 40 elements to obtain a gain of 16 
dB as a convenient number. To steer the beam, we change the phase of each element in 
accordance with Equation (13). 

Figure 17a shows the realization of a random array of 40 half wave dipole 
elements distributed over a 2.5Ax 2.5A area, and in Figure 17b the coordinate system is 
presented for reference. The minimum spacing between elements is set to O.lX, but would 
need to be increased to realize this example with the MICAz mote. The O.lX spacing 
found in some parts of this realization is much denser than a traditional planar array of 
the same area. As before, R{0,(t)) is used to visualize the coordinate system by 
representing an arbitrary angle. The cyan circle depicts the locus of observation points, 
beginning at 6 = 90° and with <p = 0° on the x-axis. Figure 18 shows the gain pattern of 
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the array steered to 0o=9O°, 0o = 0°. Figure 19 shows polar plots of the gain patterns for 
the first six beams to envision the scanning performed by the array. 
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Figure 17. Searching WSN: 

(a) Realization of a 40 Element, 2.5 A x 2.5 A Random Planar Array 

(b) Coordinate System with Scan Plane at 0 = 90° Shown in Cyan 


In Figure 18, we see the main lobe with a peak gain of 16.7 dB, comparable to the 
theoretical gain of 16 dB for a 40 element array. A random shape characterizes the side 
lobes, with the expected pedestal nature shown in Figure 12. The average side lobe level 
for all scan angles is at 17.1 dB below the main lobe. The half power beam width in 
Figure 18 is 22°, in close agreement with the expected HPBW of 50.8°/2.5 = 20°. 


42 






















Figure 18. FI-plane Gain Pattern for a 40 Element Random Planar Diploe Array 

with Beam Steering Angle 0o=9O°, (po = 0° 


Figures 19a-f show the gain pattern for six conseeutive beams in a polar plot 
format. Figure 19a is the same data as Figure 18 repeated in a polar format. 
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Gain Pattern, Scan Angle = 0.0 



Gain Pattern, Scan Angle = 80.0 
105 - 90 ' 75 - 



Gain Pattern, Scan Angle = 40.0 



Gain Pattern, Scan Angle = 120.0 
105 - 90 ' 75 - 



Gain Pattern, Scan Angle = 160.0 
105 - 90 - 75 - 



Gain Pattern, Scan Angle = 200.0 
105 - 90 - 75 - 



(e) 

Figure 19. H-plane Patterns for a Random Planar Array of 40 Dipoles as Beam is 

Scanned From 0 to 200° in 40° Steps 
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We see that as the beam is steered in each direction, it retains approximately the 
same gain and average side lobe levels on each plot. 

It should be noted that in this example we used a 20° beamwidth to reduce the 
number of beams required to scan the entire horizon, and thus the total power to complete 
one sweep. However, to obtain this beamwidth, the array size is limited to 2.5A x 2.5A, 
which in turn limits the total number of elements that can physically fit within that area. 
A reduced number of elements results in a lower peak gain. Thus a tradeoff is required 
between the total number of elements (peak gain) and the beamwidth (size of the array 
and total energy consumed to scan the horizon). 

These plots demonstrate that the beam of a WSN composed of uniformly 
distributed, randomly placed sensor nodes can be used to form a beam and rotate that 
beam across the horizon. The characteristics of the antenna beam are consistent with the 
theory presented in this chapter on random antenna arrays. This completes the discussion 
on the transmitting or searching WSN. As the beam from the transmitting array scans the 
horizon, it is assumed to illuminate a receiving WSN, which must recognize and reply to 
this signal. 

In this chapter, we presented results demonstrating control over the half power 
beam width, and gain and direction of the radiation pattern for a random planar array. 
We demonstrated that we could select a wide beam for the search pattern and 
simultaneously select the gain. This beam was steered in a lighthouse fashion simulating 
the search for and hopefully illumination of an adjacent WSN. 

We now turn our attention to the receiving array. 
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IV. DIRECTION OF ARRIVAL 


The approach of this chapter is essentially from the viewpoint of the receiving 
array. When the receiving WSN detects a signal from the searching WSN, it needs to 
know the relative geometries of the two WSNS in order to form its array and then send a 
reply. The location information of the sending WSN is contained within the message 
signal transmitted, but the direction of arrival (DOA) must be calculated. This chapter 
establishes the need to blindly determine the DOA of the received signal and presents a 
short survey of techniques available in the literature. A spread spectrum space division 
multiple access technique is introduced and demonstrated first with a linear array, and 
then applied to the random planar array. A simulation of the message extraction from a 
received signal is presented. The strengths of this technique as well as the importance of 
the spread spectrum pseudorandom noise codes are explained. 

A, DETERMINATION OF THE RELATIVE GEOMETRIES 

Here we consider two of the WSNs depicted in Figure 3, WSNs “A” and “B.” We 
consider that WSN A has sensed a phenomenon and needs to pass this data to the next 
WSN, say WSN B. In Figure 20, WSN B has been illuminated by the searching WSN A 
and now needs to determine the DOA of the signal. The receiving WSN B will form a 
random planar array using its sensor nodes to analyze the incoming signal and determine 


the DOA. 




WSN A 


WSNB 


Figure 20. Signal From WSN A Illuminating WSN B 
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The receiving WSN will then use this calculated DO A and the location 
information encoded within the transmitted signal to determine the relative geometries 
and form the optimal return beam to the searching WSN. 

1. Link Initiating Signal 

In Figure 20, the receiving array has no prior knowledge of the searching WSN 
from which the initiating signal is coming. It will have to determine the direction of 
arrival by processing the signals received at each of the sensor nodes. Section B of this 
chapter discusses this in detail. 

Some sensor nodes have the capability to operate in a power saving or “sleeping” 
mode. This thesis does not address the sleeping or transition of modes for each of the 
sensor nodes. Instead, we assume that either the search beacon contains a “wake up” 
signal or some other mechanism has allowed all sensors in the receiving array to be 
turned on for receiving signals. 

2. Return Signal 

We assume that in addition to the location of the sending array, the message in the 
search signal contains the “name” or “address” of the searching WSN, e.g., its MAC 
address and its lat/long coordinates. The receiving WSN is then in a position to solve its 
half of the link equation and form a return signal. For example it may have determined 
that WSN A at bearing 040° and range 500 m has requested a link. When it forms its 
reply or return message to the originator, it must include its lat/long coordinates as part of 
the initial handshake. The first task, however, is to determine the DOA of the query. 


B, DOA DETERMINATION 

Godara [10] provides a summary and comparison of various direction of arrival 
methods including spectral estimation, minimum variance distortionless response 
estimator, maximum likelihood, maximum entropy, least mean squares, conjugate 
gradient. Multiple Signal Classification (MUSIC) and Estimation of Signal Parameters 
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via Rotational Invariance Technique (ESPRIT). MUSIC and ESPRIT are also known as 
eigenstructure or subspace methods and can be computationally intensive [10]. ESPRIT 
and most versions of MUSIC require periodic arrays [20], thus making them unsuitable 
for the WSN in our problem. Time difference of arrival (TDOA) is another technique 
used in many applications [9]. The direction of arrival problem is also known as angle of 
arrival, spectral estimation, or bearing estimation. 

1, The Requirement for a Blind DO A Algorithm 

Prom the perspective of the receiving WSN, when a signal is received from the 
searching WSN, its DOA is unknown. The receiving WSN must, therefore, determine 
the direction from which the searching WSN is transmitting without the aid of a reference 
signal as is used by some DOA methods. DOA determination methods that do not use 
prior knowledge of signal direction or some sort of reference signal to form an initial 
estimate for the DOA are termed blind DOA methods [10]. 

A new method based upon spread spectrum techniques for space division multiple 
access (SDMA) applications that does not involve iterative matrix solutions, reference 
beams or specific array geometry has been presented and described in a patent by Elam 
[25]. Because this method does not possess the computational intensity of the above 
methods, is suitable for random arrays and has the capability for blind DOA 
determination, we chose this method for the receiving array. The rest of this chapter 
discusses and applies the SDMA technique to our problem. 


C. SPACE DIVISION MULTIPLE ACCESS RECEIVER 

The direct sequence spread spectrum (DSSS) technique proposed by Elam [25] 
and as described by Gross [20] is summarized in this section. Analogous to the code 
division multiple access (CDMA) receiver, this space division multiple access receiver 
(SDMA) spatially separates received signals using DSSS techniques. The method 
involves chipping the phase of the received signals at each array element with individual 
spreading sequences. In computer memory identically chipped virtual signals are 
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generated from a virtual reeeiver array for / expeeted direetions of arrival. The reeeived 
signals are then summed and a quadrature eorrelation is performed between the reeeived 
signal and eaeh of the / expected DOA virtual signals. Correlation values (i?/) exceeding 
a threshold are identified as signals received from an expected DOAi. From the phase of 
the correlation Ri, the message signal information may be extracted. Interfering signals 
are not well correlated and thus minimized [20]. This is similar in concept to a matched 
filter bank in a code division multiplexing receiver. 

The time sample of the received signals is of length T. With this time sample, M 
pseudorandom chips are applied to the received signal at each receive element. Each 
receive element has a separate pseudorandom noise (pn) spreading code. The spreading 
sequences can be thought of as pseudorandom array weights. If viewed in the context of 
anteima array weights, during the time period T, the array has been “steered” by M sets of 
pseudo-random weights. Over all M chips, the average pattern approaches that of an 
omni-directional array [20]. Viewed another way, the “look direction” beams are created 
through the correlation with the virtual beams. Thus, steering the “look” direction of the 
physical array is done through the software that creates the virtual array expected DOAs. 


We now discuss the receiver in greater detail. Figure 21 illustrates the receiver 
(after [20,25]). In the following sections we examine the three main components of the 
receiver; the actual array, the virtual array and the correlation and phase detection. We 
then describe the method for extracting the message in the received signal. 
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SDMA Receiver (after Gross and Elam) 
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Figure 21. 















1 . 


The Actual Array 


The actual array represents the physical receive antenna elements in the receiving 
WSN with a given geometry and L signals arriving at the array. The signals are assumed 
to have been phase shift keyed (PSK). The L phase modulated incident signals .s, (t) 

arrive at the receiver and form the received signal Xn(t) at each of the N antenna elements. 
The N elements need not be uniformly spaced, but can be arbitrarily located. Indeed, 
randomly placed elements, in a two or three-dimensional array helps the performance by 
producing independent phasing at each element for each expected angle of arrival. Each 
received signal x„(t) is then spread or chipped with a pseudorandom noise (pn) code 


associated with that antenna element, producing the spread signal (t), and finally the 
summed signal r(t). The signal is sampled over a time period of t to t+T. The received 
signal at each array element can be expressed as: 

L 

= ( 20 ) 

;=i 

where for a uniform linear array with k as the wavenumber and d the distance between 
elements, the array factor contribution , for element n from signal / is given by: 


^ _ jkd{n-l)sinei 

^nj ~ ^ 

and the baseband phase modulated received signal of interest Sf (t) is given by: 


( 21 ) 


s,{t)= (22) 

From Equation (20), we can express the received signal at each of the receive elements as 
an A X 1 column vector: 


x^it) ... x^(or 

For N elements and L signals, the received signal can be expressed in matrix form as: 

x{t)-As{t) (23) 

where A is the NxL matrix of the steering vector elements, and s{t) is an Axl 
column vector of the incident signals Si(t). The matrix A can be compactly expressed as 


a 1X L row vector of the steering vectors (0j): 
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x{t) = \a[{e;) a[{e^) ■■■ a[{Gj\ 


(24) 


'/.(O 

g?2(0 


or A can be expressed showing the expanded steering veetors in order to see the 
individual contribution of the ri^ element and signal; 



1 

1 

1 

1 

1 _ 

II 

^jkdsinO^ 

^jkd sin Q 2 

^ jkd nine L 

g?2(') 




^jkd{N-l)smd^ 

1- 

'TX 

1 _ 


For example, for L = 1, x{t) is given as: 


v(0 = 


gtfi (0 

gtC’iCOgjWsine, 


The pn chipped signal r{t) is described in matrix form as: 


(25) 


(26) 


where 


r{t) =w(t)x(t) 


(27) 


>v(o = K(o wj(o ••• 


(28) 


is the 1X row vector of the spreading codes. These spreading codes Wn(t) will be used 
to phase modulate the array outputs Xa(t). Each vector Wn(t) contains M chips of length x, 
chosen such that the length of the vector is Mt - T , the sample time. The number of 
chips M should be chosen such that it is greater than the Nyquist rate of the received 
baseband signal modulation [20]. This implies the phase modulation of the received 
baseband signal is nearly eonstant over T. The time-bandwidth product of the sample 
time T and the message signal bandwidth B should be TB < 0.25 [20]. 


In spread spectrum applications, there are many types of pseudorandom noise 
spreading codes available including maximal-length sequences or m-sequences. Gold 
codes [26], and Golay eodes [27]. Walsh-Hadamard codes [27] and Welti codes [28] can 
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be used as orthogonal eovering spreading eodes. In our implementation we have ehosen 
Walsh-Hadamard eodes beeause of the simplieity of implementation. 

The eross eorrelation of eodes within the same family provides one measure of 
how well a spreading code works [12], Codes that have a minimum cross-correlation 
will cause a minimum interference with other codes because signals spread with one code 
will not correlate with another code and be suppressed [26]. The orthogonality of the 
Walsh-Hadamard or Welti codes provides separation for each of the channels, here 
represented by receive elements which separate out the signals from each direction [20]. 
Welti codes have been reported to produce the best results in this application [29]. 

Figure 22 shows the effect of the spreading codes on the received signal at one of 
the antenna elements. Here we show the received signal phase one of the antenna 

elements in blue, the Walsh-Hadamard spreading code Wn in red, and the real part of the 
chipped signal r^(t) in purple. The chipping process is repeated at each of the N receive 
elements, each element having its own individual Walsh-Hadamard code. 
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Figure 22. Received Signal Phase Walsh-Hadamard Spreading Code Wn, 

and the Real Part of the Chipped Signal r„(t) 


Key to this technique is the formation of the virtual array of signals that “look” in 
the expected directions of arrival. The generation of the virtual signals is similar to that of 
the actual array. We continue with the description of the SDMA receiver by examining 
the virtual array component. 

2, The Virtual Array 

Figure 23 graphically shows this component of the receiver. Within computer 
memory, a “virtual array” is configured according to the geometry of the physical array, 
containing the same number of elements. The virtual array is steered to an expected DOA 
by computing the array factor for that expected direction 0,. The expected direction Oj 
can also be thought of as the beam steering angle Oo . The array factor for each of the N 
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elements is chipped with the identical spreading codes as the physical array. This is 
repeated and stored for each of the / expected directions of arrival 0,. The / expected 
directions can be chosen to avoid known interfering directions, prior knowledge of the 
incoming signal, multipath [20,25], etc. The expected directions cannot be chosen any 
finer than physically allowed by the actual array. 


Virtual Array 




Figure 23. Expected Signal Generation by the Virtual Array in Memory 


For each of the / expected directions, we can express the summed virtual signal 
as: 


v,(t)= w(Ofi;(0,) (29) 

where the steering vector for each expected direction is given by the V x 1 column 
vector: 



1 

jkdsixi 6: 

e 


and the spreading codes are the same as before: 

w{t)^\w^{t) w^it) ••• 


(30) 


Having formed both the received signal and the array of virtual signals, the 
receiver now performs a correlation computation to establish the DOA of signals incident 
to the receiver. 
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3, Signal and Phase Correlator 

The signal and phase correlator correlates the actual received signal r(t) with 
each of the / virtual signals v,. (t), over the sample time T. Figure 24 illustrates the 
operation. 


r{t) 



Figure 24. Signal and Phase Correlator 

The correlation computation is expressed as: 

= r{t)v*{t)dt (31) 

where Ri is the complex correlation value for direction 0 ^. The magnitude of Ri will be 
plotted for each angle 6, becoming the array factor for the signals present. The threshold 
detector depicted compares the magnitude of Ri against a predetermined value and 
declares a signal present or not at this angle. If a signal is present, the data is extracted as 
explained in the next subsection. This is similar to a matched filter or correlation detector 
operation [27]. 

The correlation values, when the pn codes for the actual and virtual signals match, 
become the array factor for the antenna array in the direction of the virtual signal. To see 
how this works, we examine the correlation in more detail. From Equations (22) and (27), 
the actual signal received at the antenna array r(t) with N array elements is expressed as: 
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( 32 ) 


^(0=X 


^i(n-\)kd sinSo 


n=\ 


where the spreading codes Wn are modulating the phase of the received signal. The signal 
generated in memory for the expected direction and modulated with the identical 
spreading codes is from Equation (29): 


V; 


N .K 


(33) 


Substituting Equations (32) and (33) into (31) and noting that the signal has been sampled 
in time through the spreading process gives the expansion: 


R 


A' 

-IS 


m=l 


_ .n 


(34) 


Rearranging the element summation yields 


N N 




^j(n-l)kd sin 2 sine,. ^ i mi*) 


(35) 


f n=\ m=\ 


Since the array factor terms are independent of time, they can be moved outside 
the integral. Recall that we chose the time sample period T such that the phase ^(t) of the 
signal s(t) was approximately constant over the sample period and can be approximated 
as its average ^ ; thus, it can be moved outside the integral as well: 


«=1 m=l t 


5{n-m)=- 


1 if n=m 
0 if n^m 


(36) 


The integral yields a delta function 5{n-m) because and are orthogonal [12]. 
Thus Ri reduces to: 


_ ^gy(«-l)<^(sin0„-sine,) 


(37) 


Equation (37) is equivalent to the array factor, here multiplied by the constant 
phase of the received signal. Consequently, the correlation can be expressed in terms of 
magnitude and phase as: 

R,=\R\e^P' (38) 
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where |/?, | is the magnitude of the eomplex valued correlation and p, expresses the phase 

of the received signal. The magnitude of the correlation contains the DOA information 
and the phase term contains the message information. 

4. Message Extraction 

The phase of the signal s(t) can also be estimated, and the message data extracted 
using this SDMA method. Over the time sample period T, the phase of signal s(t) was 
assumed to be approximately constant. Thus, when a signal is detected from direction z, 
the phase and thus the message can be approximated as p,. The correlation integral 
averages the received signal phase over T [20], therefore 

— 1 ft+T 

= um (39) 

where ^ is the average signal modulation during the time period T at 6i. By calculating 
p, at each angle where a signal has been determined present, we can simultaneously 

extract more than one message signal from the correlation data. The process is described 
in further detail in [20, 25]. 

5, The Importance of the Spreading Codes 

One may ask what advantage do the spreading codes provide to the overall 
receiver? The answer, in short, is that the weight functions or the chipping sequences 
provide orthogonality of each received signal to all of the incorrect DOA virtual signals. 
To examine this, we consider an example with A = 4 elements; the signal is incident from 
direction Oinc and Wn, k and d are defined as before. The received signal is then: 

r{t) = ] 

Recalling that the time sample period T was chosen such that s(t) was 
approximately constant in phase, the leading exponential can be suppressed as a constant. 
If r(t) is not chipped, the time sample is reduced to a single sample point in time, 
and the expression for r(t) becomes the summation: 

N N 

= = (40) 

«=1 n=l 
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where represents an array elements’ eontribution to the array factor. 

For each of the / expected directions, the virtual signal is expressed as; 

(41) 

n=\ n=\ 

where dm represents the virtual array elements’ contribution to the array factor. 

Using Equations (40) and (41) in the correlation computation yields: 

= r(Ov:(o=m 

n = l m=l tt = l m = l 

This expression is the correlation sum that results when the spreading codes are not used. 
Now substituting: 

= Msin0;„^ and \i/. = kdsinO^ (43) 

into Equation (45) and expanding the sums, the correlation sum takes the form; 

The correlation magnitude should be at a maximum when the virtual angle 0, matches 
the incident signal angle i.e., y^i„c = Wi - ^dsmd.^^= a constant. After some 
algebraic manipulation. Equation (44) can be written as; 

3 1 

^max “ ^ 2 ^inc ) + cos(2Msin) + - cos(3Msin) (45) 

While our intent was that the correlation reached its maximum when the expected 
angle matched the actual incident angle, this clearly is not the case. In Equation (45), the 
cosine and sine terms will influence the value of Ri depending upon the particular value 
of 0 for the incident signal. Thus evaluating the correlation without chipping codes will 
not give us the clear solution we are seeking. While the correlation between r{t) and 
v.{tg) should be large when Oj matches the actual received direction of arrival 0,„c, there 
is no guarantee that for any other arbitrary angle 0/, the correlation would not be large. 

More generally, we can clearly see the advantage of the chipped sequence by 
comparing the correlation with and without the spreading sequence. In the non-spread 
approach, the correlation sum contains cross terms as seen in Equation (42) while the 
SDMA approach produces the array factor in which the cross terms are zero as given in 
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Equation (36). It is clear that the orthogonal spreading codes used in the SDMA 
technique are integral to the method. 

6. Advantages of the SDMA Technique 

The strength of this technique comes from the ability to receive several signals 
and “look” in almost all directions nearly simultaneously. Because the chipping of the 
phase of the received signals and the phase offset at each of the antenna elements, this 
essentially makes the set of possible antenna patterns independent during the sampling 
period. The orthogonality of spreading codes between antenna elements allows the actual 
signals to be extracted with the eorrect direction of arrival. The signal message is 
extraeted from the phase of eorrelation sequence. Since the set of virtual signals is 
computed once and stored in memory, the computational cost is significantly reduced as 
compared to the traditional DOA methods, sueh as conjugate gradient and other iterative 
matrix solution methods. Additionally, the teehnique provides a processing gain of 

[ 25 ]. 

Having described the SDMA method we now demonstrate through simulation the 
effectiveness of the method as applied to a uniform linear array and a random planar 
array of antenna elements. 


D, MATLAB SIMULATION RESULTS 

In this seetion we demonstrate the results of the SDMA technique using 
MATLAB. First we apply the technique to an 11 element ULA and demonstrate the 
importance of ehoosing the optimum spreading codes. Having seleeted the Walsh- 
Hadamard codes we then show the affect of the number of chips used to spread the 
incoming signal. We then apply the teehnique to the random planar array and simulate 
solving for a DOA and then extracting a phase encoded message from the detected signal. 
This is a simulation of the WSN B of Figure 20 receiving a message from WSN A and 
then deeoding the location message. In the message extraction we simulate a white 
Gaussian noise channel. 
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1 . 


Uniform Linear Array 


In this simulation we use the SDMA method to ealculate the DOA of two 
simultaneous test signals. The ULA contains 11 z-direction oriented dipole elements 
spaced one-half wavelength apart and is used to detect test signals arriving from two 
different directions -30° and +45°. The array factor for this ULA is given simply as the 
sum of the array factors from each incident angle. For an N element ULA, the array factor 
for a single incident wave from direction Omc is [18]; 

sin [Nkd (sin0 - sin 0.^^ 

N sin (M (sin 0 - sin 0;„^)) 

and the complete array factor for two incident signals is the sum of both array factors; 

■W»,(«) = Af^(8) + AF^y(e) 

Again, the incident “look” angles can be thought of as the beam steering angles. Here, 
the incident beams are restricted to the xz plane. The geometry used in this simulation is 
shown in Figure 25. 




Figure 25. Geometry for ULA with Signals Incident from -30° and +45° 

The choice of spreading codes has an influence on the quality of the DOA 
determination. In Figure 26, we compare the result of the correlation using two different 
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types of spreading eodes, the Walsh-Hadamard eodes and a pseudorandom noise code 
sequence. The correlation magnitudes and the array factor are plotted against the angle of 
incidence 0. 



Figure 26. SDMA - Comparison of PN and Walsh Spreading Codes Using 32 Chips 

11 Element ULA with Signals Incident from -30° and +45°. 


The green line is the array factor. In this comparison, we used 32 chips for both 
the pn sequence and the Walsh-Hadamard codes and 4 samples of the signal per chip. A 
low number of chips was chosen in order to make the difference in spreading code choice 
more obvious. The blue line is the correlation magnitude calculated using the Walsh- 
Hadamard spreading codes and the magenta line is the correlation magnitude calculated 
using the pn sequence. In this figure, Walsh-Hadamard codes clearly work better than a 
pn sequence. The pn codes are adequate at determining the DOA of the signal, but the 
sidelobes do not match closely with those of the theoretical array factor. For a random 
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array, where the sidelobes do not decrease as strongly as the ULA, this suggests that the 
Walsh-Hadamard codes are a better choice. 

We now examine the effect of the number of chips used. In Figure 27, we 
increased the number of chips used from 32 to 64 with the same number of samples per 
chip, array geometry and incident signals as in Figure 25. The array factor for the ULA 
is plotted in green and the magnitude of the correlation 7?, is plotted in blue. We see from 
this result that the SDMA technique produces a correlation in excellent agreement with 
the theoretically produced array factor, clearly showing the DOA of the two incident 
signals at 0 = -30° and +45° and accurately matching the mainlobe as well as all the 
sidelobes. 



Figure 27. DOA Determination of Two Incident Signals Using SDMA on an 11 Element 

Array with One Half Spacing and 64 Chip Walsh-Hadamard Code 

From Figures 26 and 27, we see that the choice of spreading code and the number 
of chips are important. In comparison to Figure 26, we see that 64 chips produce a much 
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closer match to the theoretical array factor than 32 chips when using Walsh-Hadamard 
codes on the same array. Additional comparisons (not shown) reveal that increasing the 
number of chips to 128 or greater had only marginal increase in the correlation magnitude 
as compared to theory. It appears that 64 chips are adequate for DOA determination. 

Having demonstrated that the SDMA receiver performance agrees with theory for 
the uniformly spaced linear array, we now extend the technique to the problem of the 
planar array of randomly placed elements. 


2, MATLAB Simulation Results - SDMA and Random Planar Array 


In section A of this chapter, we discussed the need for the receiving array to 
determine both the DOA of the searching signal and the range to the searching WSN. 
Here we apply the SDMA technique to the random planar array of dipoles to determine 
the DOA of a test signal and compare the gain pattern from the SDMA technique to that 
of the array as calculated in Chapter III. Having found the DOA, we simulate the 
decoding of a sample message representing bits of a data stream from the searching 
WSN. The method for calculating the contribution of each array element to the 
correlation value for a random planar array is exactly analogous to that of Chapter III; 
instead of using a geometric expression for the location of a receive element we use the 
exact location as was done in Equation (14). We modify Equation (36) for the random 
planar array as follows: 


R. = 

n=l m=l 




■ 2 e 2 dt 


(46) 


where the element phase y/^( 0 o, 0 (,) is now determined by its x-y coordinate and the 


steering angle ( 0 o’ 0 o) given by: 


V^« (^0 ’ 00 ) = (sin 0 cos (^ - sin 00 cos ^ 0 ) + (y„ sin 0 cos 0 - sin 00 cos 00 ) 


a, DOA Determination 

We now apply the SDMA technique to a planar array of randomly placed 
sensors. Eigure 28a depicts an overhead view of the elements of a 5A x 5A sensor WSN 
containing 30 elements, randomly distributed with the array scanning the x-y plane. 
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Figure 28b depicts the simulated test signal incident from 6 = 90° and (j) = 30°. In this 
diagram is shown to establish the coordinate system and indicates the direction 

of incidence of the test signal. 




He 

- » - 

H' 


^ , 

H^ 

.. 

if: 

• 

He 

He 

He 

ii 

H 

H 

.#.. 

He 

— 


(5 

- ^ 

He ^ 


X* wavelengths 



Figure 28. Thirty Element Random Planar Array Realization and Coordinate System 


With the geometry shown in Figure 28, the gain pattern of the random 
planar array was calculated using Equations (3) and (14) and compared to the SDMA 
gain pattern calculated using Equations (3) and (46). Equation (3) is the element pattern 
of the half wave dipole. In Eigure 29, the green line shows the gain pattern of the random 
planar array with the incident test signal; it is the signal the SDMA technique must locate. 
The SDMA gain pattern is presented in blue and shows excellent agreement with the 
calculated gain pattern. 
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Figure 29. Gain Pattern of 30 Element Random Planar Array, Steered to 30° with SDMA 

Estimation of Incident Signal Pattern. 

As shown in Chapter III, the random array produces a main lobe in the 
steered direction with the sidelobes displaying the characteristic random nature below a 
pedestal. The 30 elements produce a peak gain of 17.8 dB, the beamwidth is 8 degrees 
and the highest side lobe rises to 8.3 dB below the peak. These values are in agreement 
with the theoretical results presented in Chapter III. The SDMA pattern accurately shows 
the main lobe at 30° and all of the sidelobe peaks. The angle where the nulls of the gain 
pattern occur are accurately shown with the SDMA technique, but the absolute values do 
not exactly agree with the theoretical values. 

These results validate the SDMA technique as a method for blindly 
determining the DOA when applied to the random planar array. It shows excellent 
agreement between the theoretical array patterns and those produced by MATLAB. This 
will allow the receiving WSN to almost instantaneously determine the DOA of the 
searching WSN. 

Having determined the DOA of the searching WSN, the task of the 
receiving WSN is to extract the location information from the message signal sent by the 
searching WSN. 
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b. Message Extraction 

One of the primary goals of the reeeiving array is to extraet from the 
transmitted signal a message containing the location of the searching WSN. In order to 
demonstrate this, we will simulate the baseband message signal as a bit sequence and 
include it in the signal incident from 30°. The goal is to transmit an 8 bit sequence (Oil 
0 10 11) mapped to (1 -1 -11-11 -1 -1) as the sample message. The signal phase ^(t) 
was generated with additive white Gaussian noise (AWGN) with a 20 dB signal to noise 
ratio to simulate channel interference. Here MATLAB measures the power in l^(t) and 
adds AWGN at 20 dB. The signal to noise ratio is the ratio of the power in the signal ^ (t) 
to the power in the AWGN. 

The results of this simulation are shown in Figure 30. We see in blue, the 
instantaneous baseband signal l^(t) with AWGN. In green is the phase term of the 

correlation sequence Ri at the incidence angle 0 = 30°. The time sample of the received 
signals is of length T and shown in red. For simplicity, we send one bit per time period T, 
although this need not be the case [20,25]. We see that the signal bit sequence can be 
accurately extracted from the recovered estimate of the phase. 
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Figure 30. SDMA Message Extraetion of Sequenee [01101011] 
from Noisy Baseband Signal, SNR = 20 dB 

In order to test the strength of the SDMA teehnique, we increased the noise level 
such that the SNR was now -5 dB. Figure 30 shows the received baseband test signal in 
blue encoding the same bit sequence at the higher noise level. The recovered signal is 
shown in green. Although the noise level appears very high, the SDMA technique is still 
able to extract the data in the incident signal. This suggests that the receiving WSN is 
capable of extracting the message bits under high noise conditions. 
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Figure 31. SDMA Message Extraetion of Sequenee [01101011] 

from Noisy Baseband Signal, SNR = -5 Db 


The direetion of arrival for a signal from a searehing WSN is determined by the 
magnitude of the eorrelation sequenee and the message itself is eontained within the 
phase of the correlation. The bit stream is the message containing the “name” and 
“location” of the WSN searching for the adjacent WSN. Thus, the two key pieces of 
information required to establish the link, the bearing and range can now be determined. 
The latitude and longitude location information can be converted into a range using 
Equations (18) and (19). 

The results shown here are preliminary to demonstrate the effectiveness of the 
proposed scheme. Eurther in depth simulation and comparison to gain a full 
understanding and scope of the proposed technique is required. 

In this chapter, we established the need for the receiving WSN to blindly 
determine the DOA of an arriving signal. A new SDMA receiver was described and its 
effectiveness is demonstrated through the simulation of a random planar array. We 
showed that the method works better for Walsh-Hadamard spreading codes than for pn 
codes generated by a pseudorandom number generator and that a 64 chip sequence is 
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sufficient to accurately establish the DOA and extract the message. Using a random 
planar array, we validated the technique by determining the DOA of an incident test 
signal and then extracting a message data stream from this signal in a noisy channel at 
SNRs of 20 db and -5 dB. From the DOA and the location information within the 
arriving signal, the WSN has sufficient knowledge to respond to the query of the 
searching WSN. 

Once the receiving array has knowledge it is being queried for a link, the task is 
then to focus its array towards the searching array, calculating the necessary parameters 
required for the link. This is done via the link budget analysis as was described in 
Chapter II, thus completing the link between two WSNs. 
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V. CONCLUSIONS 


In this chapter, we present a summary of the work done, highlight the signifieant 
results and suggest topies for further researeh. 


A, SUMMARY OF RESEARCH 

In this thesis we explored methods for forming an Over The Horizon (OTH) 
eommunieations link of Wireless Sensor Networks (WSNs) by enabling eaeh WSN to aet 
as a smart antenna array. Methods for establishing the OTH eommunieations link via 
beam forming and direet sequenee spread speetrum Spaee Division Multiple Aeeess 
(SDMA) were presented and modeled using MATLAB. Methods for forming a seareh 
beam for a transmitting WSN and determining the direetion of arrival of the reeeived 
beam for a reeeiving WSN were presented. 

We examined the searehing WSN in whieh the individual sensor nodes operate in 
eoneert as a smart antenna array. The WSN then forms and transmits a seareh beam 
around the horizon in order to illuminate an adjaeent WSN and establish a link. Each 
WSN was modeled as a uniformly distributed random array of sensors nodes. Gain and 
beamwidth of the random planar arrays were caleulated and validated using MATLAB 
simulation and eompared to theoretieal results. Through MATLAB simulation, the half 
power beamwidth (HPBW) of a random planar array as a funetion of WSN size and the 
gain as a function of the number of elements for a fixed size array were studied. A 40 
element random planar array, eomparable to a 36 element periodie planar array of similar 
size, was simulated to demonstrate the search capability. . 

We examined the proposed eommunieations link from the perspeetive of the 

reeeiving array. This array needs to determine the direetion of arrival (DOA) and 

loeation of the searehing WSN in order to form an optimum link without a referenee or 

pilot signal. A SDMA reeeiver eapable of determining the DOA and loeation of the 

searehing WSN was deseribed and its effeetiveness demonstrated through the simulation 

of a random planar array. Using a random planar array, we validated the SDMA method 
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by determining the DOA of an ineident test signal and then extraeting a message data 
stream from this signal in a noisy ehannel. From the DOA and the loeation information 
within the arriving signal, the WSN has suffieient knowledge to respond to the query of 
the searehing WSN. 


B, SIGNIFICANT RESULTS 

We investigated the random planar array performanee in terms of the half power 
beam width, the gain of an array of N elements and side lobe performanee. The objective 
was to apply these results to wireless sensor networks to first form a search beam to 
locate adjacent WSNs, then form a narrow high gain beam to form a communications 
link. The particular WSNs are assumed to be fixed in location, but the inter-element 
spacing follows a uniform random distribution within a defined boundary. We 
demonstrated that we could select a wide beam for the search pattern and simultaneously 
select the gain. To achieve a given beamwidth, we can choose a physical size, likely a 
sub-section of the WSN under consideration. Within this physical area, we have control 
over the gain by choosing the number of elements employed in the array. Alternately, for 
a given required gain level we can thin the array (or reduce density) without significant 
loss of gain or the effects of grating lobes. It does not matter which particular sensor 
nodes within the WSN are used as long as, on average, they are randomly distributed 
within the physical area chosen. This is important if the objective is to distribute the 
burden of transmitting, and thus energy consumption among the nodes within the WSN. 

We established the need for the receiving WSN to blindly determine the direction 
of arrival (DOA) of an arriving signal. From the DOA and the range information 
encoded within the arriving signal, the WSN has sufficient knowledge to build its reply to 
the query of the searching WSN. Using a random planar array, we utilized a spread 
spectrum base SDMA technique to determine the DOA of an incident test signal and 
simulated extracting a message data stream from this signal in a noisy channel with SNRs 
of 20 db and -5 dB, accurately extracting the data. 


72 



C. TOPICS FOR FURTHER RESEARCH 

In the process of carrying out this work many issues for further research came to 
light. This section highlights some of these topics. 

In this work, we demonstrated a SDMA method to determine the DOA of a signal 
in two dimensions assuming vertical polarization. Indeed, the SDMA method shows 
promise for much greater capabilities. In this thesis, the virtual signals were only used to 
determine the DOA across the horizon but could be extended to include the hemisphere. 
In this manner, the graphically depicted UAV of Figure 4 could quickly be located, or in 
the reciprocal instance, the UAV could quickly find WSNs. Once located, the look 
directions or virtual signals used for tracking quickly collapses to just a few vectors. This 
type of three dimensional pattern covering the hemisphere above the array is often 
referred to in the literature as a “pin cushion” pattern. The SDMA technique can also 
search for multiple or elliptical polarizations by creating and summing the orthogonal 
polarization virtual signals [25]. 

In this work, the spreading was applied only to the receiving WSN to determine 
the DOA and extract the message. Alternately, the searching WSN could spread its 
transmitted search signal in the manner described in Chapter IV on SDMA for the 
receiving array leading to further processing gain. 

For the SDMA receiver, we assumed that only the direct signal reached the 
receiver. For a terrestrial communications link, there would certainly be a multipath 
signal fading [16]. Future work may investigate the receiver’s capability to reject 
multipath signals and increase the communication link range or bit rate. 

In this work, we assumed that the location of each sensor and that the location 
error could be estimated. Appendix A addresses the methods of self-localization of the 
array, however the impact of these errors to the SDMA technique were not addressed in 
this thesis. Future work may determine the impact of localization errors on the SDMA 
technique. Additionally, future work may address methods for combining the signals 
received at each of these individual sensor node locations. 
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The path loss model in this work used a wireless LAN model based on antenna 
heights of 1 m or above. A better path loss model would give a more accurate estimate of 
the margin required to form the communications link between two adjacent WSNs. 
Future work may seek to develop path loss models for ground level transmitter-receiver 
pairs based on empirical data using actual sensor motes under a variety of operating 
conditions. Applications for linking together WSNs are only limited by the imagination. 
The WSNs need not be random or fixed in place on the ground. For example, one WSN 
may be a periodic array placed on the side of a building and the other could be an array 
placed on the top or side of a truck. As the truck passes a station, the fixed station would 
search for the truck’s array, establish a communications link and pass information 
wirelessly to the fixed node. Since the link is established autonomously, the truck 
doesn’t need to stop, or indeed perform any action. A future effort may examine such 
applications that have potential military value. 
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APPENDIX A. WIRELESS SENSOR NODE LOCALIZATION 


A basic requirement of ereating a smart antenna from a WSN is the immediate 
task of forming a loeal network of these nodes, identifying eaeh other’s loeation and 
relative geometry. This map of the geometry is ereated during the self organization of the 
WSN. Onee the loeal geometry of the WSN is determined, it will then be possible to 
form a beam eolleetively and seareh for other WSNs. This seetion deals with 
determining the loeal wireless sensor network geometry. 

It is assumed that the sensor nodes are deployed or dispersed in a random fashion 
in the area of interest and that all are properly oriented, i.e. they are not laying on their 
side or upside down, but have the transmitting antenna pointing vertieally. For the 
purposes of this work, the sensors in a network are assumed to be randomly distributed 
with a uniform probability density funetion (pdf). 


A. LOCALIZATION METHODS 

So important is the need to determine loeal geometries of an ad hoe wireless 
network, extensive researeh has been eondueted on this topic and the eorresponding 
literature covers a multiplieity of techniques [4,30,31,32,33]. Here we present a brief 
deseription of the methods for determining the geometry of an ad hoe wireless sensor 
network. Due to the available methodologies and advaneing teehnology, we assume in 
this thesis that eaeh loeal WSN is self aware of its absolute loeation and orientation as 
well as the geometry of the individual nodes that eompose the network. 

Self organization allows the WSN to determine its own geometry. These 
teehniques typically estimate distanee and angle between nodes, eommunieate them to 
the neighboring nodes, then adjust or rotate individual eoordinate systems to obtain a 
single global eoordinate system. Magnani and Leung [4] diseuss one sueh algorithm. 

Loealization methods generally fall into 4 major categories and are deseribed 

here. 
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1 . 


GPS Localization 


All the sensor nodes have on-board GPS reeeivers and simply eommunieate this 
to adjaeent nodes, forming the table of loeations and or separation distanees. 

2. Signal Techniques 

These teehniques rely on nodes estimating the loeation of other nodes by 
measuring a signal parameter [30]. One sueh method measures the reeeived signal 
strength as eompared to a known transmitted value ealled Reeeived Signal Strength 
Indieator (RSSI) or Relative Signal Strength (RSS). Other teehniques are Time of Arrival 
(ToA), Time Differenee of Arrival (TDoA) and Angle of Arrival (AoA) estimation. For 
the AoA ease, the sensor node must be equipped with an antenna that is able to measure 
direetion. 


3, Anchor Nodes 

Anchor nodes or “base stations” of known position are distributed across a sensor 
network [30,31]. Other nodes’ positions are triangulated based on the known locations of 
the anchors. Careful selection of the base stations’ position distribution and precise 
knowledge of their location are crucial to minimizing position error propagation as 
subsequent node locations are calculated [31]. Additionally, the triangulation error in 
locating a particular node can be reduced by dispersing the base stations at the perimeter 
edges of the sensor field as opposed to picking 3 base stations that are close together. 

4, GPS Free Methods 

Numerous references are available for the localization of WSN in the absence of 
GPS data [4,32,33]. These methods require a signal technique and the determination of 
local and absolute coordinate systems. Reference [31] shows that multi-path errors 
dominate the measurement noise and subsequently tend to be the main cause of error in 
range estimation. 

Ji and Zha [30] discuss the method of multidimensional scaling as technique to 
mitigate the problems of cumulative measurement error. This technique builds a series of 
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local maps of adjacent sensors along a route from a starting anolior to an ending anehor. 
This method may be of use for the larger problem, diseussed in Chapter II on building the 
OTH network of WSNs. 

In this work, we assume the following about the WSNs: the WSN is loeation 
aware with respect to its absolute loeation and orientation; the relative positions of eaeh 
of the sensor nodes eomprising the network are known; the sensor nodes are uniformly 
distributed within a definable boundary; and the sensor nodes are able to synehronize 
internal cloeks either to a master node or to a GPS referenee. 
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APPENDIX B. MATLAB CODE 


This appendix contains selected MATLAB code used to generate the figures in this 
thesis. 

A, Comparison of Power Pattern of Linear Array and Average Power Pattern 
of the Same Size Random Array 

This code was used to generate Figure 12. 


% Chris Taylor Modification to Array factor code presented as a solution to 
% Homework 4.1 of below textbook. 

% This code produces both the Array factor for a uniform linear array and 
% the envelope from a random array of the same size. 

% This code also used in the Figure 9 on grating lobes. 

% Feb 08 Naval Postgraduate School, Monterey CA 

% 

% Algorithm for ULA by Frank B. Gross "Smart Antennas for Wireless 
% Communications", c 2005, Chapter 4, Homework solutions. 

% 

% Matlab m-file for Problem 4.1 

d = [.5 0.5 0.5]; % d/lambda 

%N = [888]; theta = -pi/2:.01:pi/2; 

N = [5 8 20]; theta = -pi:.01:pi; 

% -Create Normalized Array Factors- % 

AFnl = l/NCl)*sinCNCl)*pi*dCl)*sinCtheta))./sinCpi*dCl)*sinCtheta)); 

AFn2 = l/NC2)*sinCNC2)*pi*dC2)*sinCtheta))./sinCpi*dC2)*sinCtheta)); 

AFn3 = l/NC3)*sinCNC3)*pi*dC3)*sinCtheta))./sinCpi*dC3)*sinCtheta)); 

% -Linear Array Power Pattern- 

PAFnl = AFnl.*conjCAFnl); %PAFnl=10*logl0CabsCPAFnl)); 

PAFn2 = AFn2.*conjCAFn2); %PAFn2=10*logl0CabsCPAFn2)); 

PAFn3 = AFn3.*conjCAFn3); %PAFn3=10*logl0CabsCPAFn3)); 

%—Random Array Predicted Power Pattern Envelope- % 

RAFNl = CAFnl.*conjCAFnl)).*Cl-Cl/NCl)))+ 1/NCl) ; 

RAFN2 = CAFn2.*conjCAFn2)).*Cl-Cl/NC2)))+ 1/NC2) ; 

RAFN3 = CAFn3.*conjCAFn3)).*Cl-Cl/NC3)))+ 1/NC3) ; 

thetad=theta*180/pi; 

% -Plot Results- % 4-1 

subplotC311), plotCthetad,10*logl0CPAFnl), 'b' ,thetad,10*logl0CRAFNl), '-c'), gcid on 
titleCE {'N Element Linear Array Power Pattern'},... 

{' vs Average Power Pattern of N Element Random Array'}],... 

'FontSize' ,14, 'FontWeight' , 'bold') 
ylabelC '\fontsize{14}AF^2C\theta) dB ') 
xlabelC '\fontsize{14}\theta (deg) ') 

legendC'N = 5, d/\lambda = 0.5', 'Random Array Ave PP'), 
axisC[0 90 -20 0]) 

subplotC312), plotCthetad,10*logl0CPAFn2), 'r' ,thetad,10*logl0CRAFN2), '-m'), grid on 
ylabelC '\fontsize{14}AF^2C\theta) dB ') 
xlabelC '\fontsize{14}\theta (deg) ') 
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legendC'N = 8, dAlambda = 0.5', 'Random Array Ave PP'), 
axisC[0 90 -30 0]) 

subplotC313), plotCthetad,10*logl0CPAFn3), 'b' , .. . 

thetad,10*logl0CRAFN3), 'g'), grid on 
ylabelC '\fontsize{14}AF^2C\theta) dB' ) 


xlabelC '\fontsize{14}\theta (deg) ') 

legendC'N = 20, dAlambda = 0.5', 'Random Array Ave PP'), 
axisC[0 90 -40 0]) 


B, Random Planar Array Half Power Beam Width vs Array Size and 

Random Planar Array: Maximum Gain vs Array Size, Const Number of 
Elements 

This code was used to generate Figures 14 and 16 

% This code computes far-field gain pattern 
% for planar array composed of N vertical half-wave 
% dipoles at wavelength normalized coordinates [xCn),yCn)] 

% with complex ICn) currents. 

% 

% Original code for period arrays by M. A. Morgan 29 Oct 07 
% 

% Modified by C.E. Taylor Dec 07 to use Random Planar Array 
% ** This version calculates half power BW vs array size ** 

% 

% This version uses Uniformly current amplitudes CIn=l for all elem) 

% 

clear all 
clf reset 


% Setting up 3-D pattern grid 
Ntheta=91; Nphi=361; 

theta=linspaceC0,pi/2,Ntheta); phi=linspaceC-pi,pi,Nphi); 

[Phi,Theta]=meshgridCphi,theta); 

CT^cosCTheta); ST^sinCTheta); CP=cosCPhi); SP=sinCPhi); 

Rx=ST.*CP; Ry=ST.*SP; Rz=CT; % unit radial vector components 

% Half-wave dipole pattern for Iin=lA (eps handles theta=0 indeterminent) 
E0 = 60*cosCpi*CT./2)./CST+eps); 


HPBW_ne = zerosCl,50); 
Gain_ne = zerosCl,50); 
PeakSLL_ne = zerosCl,50); 
AveSLL_ne = zerosCl,50); 


% Half Power Beam Width vs array size 
% Gain vs array size 
% Peak SLL vs size 
% Ave SLL vs array size 


% -Build Random Array- 

lambda = 1.0; 
k = 2*pi/lambda; 

%array_size = 5.0*lambda; % in wavelengths (lambda) square grid 


nruns = 35; 
for run=l;nruns; 


fprintfC ' run; %2.0f\n', run); 
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for array_size = Cl;20)*lambda; 
num_elements = 20; 

Nx = num_elements; 

Ny = num_elements; 
sepn_min = 0; 

xp = array_size * randCl,num_elements); 
yp = array_size * randCl,num_elements); 

% -Stats of Random Array- 

min_spacing = 0.1*lambda; % l/10th wavelength 
while Csepn_min< min_spacing) % min spacing = .1 wavelength 

rp = Cxp;yp]'; 

sepn = pdistfrp); % dist between all elements in Random Array 
sepn_min = minfsepn); % find minimum separation 

zsepn = squareformfsepn); % convert linear sepn array to square 

[Rmin Cmin] = findCzsepn==sepn_min,l, 'first'); %returns index 
fprintfC['Min separation is %6.3f wavelengths,'... 

'between elements %2.0f and %2.0f \n'],... 
sepn_min, Rmin, Cmin) 

if Csepn_min< min_spacing) % then Move one of them 
xpCRmin)=array_size * rand; 
ypCRmin)=array_size * rand; 
end 
end 

fprintfC 'Element %2.0f Located at %8.3f, %8.3f\n',Rmin,rpCRmin,;) ) 
fprintfC 'Element %2.0f Located at %8.3f %8.3f\n',Cmin,rpCCmin,;) ) 

10=1; % assuming uniform current amplitudes 

%phid0=inputC'Enter peak gain azimuth angle in degrees from x-axis: ') 
phid0 = 0; 

phi0=pi*phid0/180; %convert to radians 

% Phasing ICNx,Ny) array elements to point array in phi0 direction 
Ip= 10 * expC-j*k*CcosCphi0)*xp+sinCphi0)*yp)); 

% Summing far-field pattern 

F=zerosCsizeCE0)); 
for n=l;num_elements 
xpRx = xpCn)*Rx; 
ypRy = ypCn)*Ry; 

F = F + IpCn)* expCj*k*CxpRx+ypRy)); 
end 

% Normalized power pattern 
F2 = absCCE0.*F)/Cnum_elements)).^2; 

DTheta=thetaC2)-thetaCl); DPhi=phiC2)-phiCl); 

IntF2=DTheta*DPhi*sumCsumCCF2).*sinCTheta))); 

Gain=4*pi*F2/IntF2; 

Gmax=maxCmaxCGain)); GmaxdB=10*logl0CGmax); 

% Horizontal pattern (the x-y plane) and 3dB beamwidth 
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n=findCphi==0); 
m=findCtheta==pi/2); 
G3dB=Gmax/2; 

GH^GainCm,;); 

mf^findCGH >= G3dB,l, 'first'); 
mp^findCGH >= G3dB,l, 'last'); 
BWH=180*CphiCmp)-phiCmf))/pi; 
GHdB=10*logl0CGH); 


% - calc Side Lobe levels - 

BW = imregionalmaxCGHdB); %MATLAB rtn that finds local maxima->BW array "l"s 
A1 = findCBW); %index pointers to local Maxima 

[bb,ix]=sortCGHdBCfindCBW)), 'descend'); % sort by highest peak->lowest 
SLL= bbCl)-bb; % Side Lobe levels, down from peak 

SLi= AlCix); % index to peaks 

SLa= 180*phiCSLi)/pi; % Angle where peaks occur 

% fprintfC'First Sidelobe Level; %6.2f below max at phi= %4.0f \n',SLLC2),SLaC2)) 

% fprintfC'Second Sidelobe Level; %6.2f below max at phi= %4.0f \n',SLLC3),SLaC3)) 

% fprintfC'Mean Sidelobe Level; %6.2f below max \n',meanCSLLC2;lengthCAl)))) 

HPBW_neCarray_size) = HPBW_neCarray_size)+BWH ; % Gain vs number elements array 

theta3db Cai^i^ay-Size) = 50.76 / array_size; % Theory 

Gain_neCarray_size)= Gain_neCarray_size) + GmaxdB; 

PeakSLL_neCarray_size)= PeakSLL_neCarray_size)+SLLC2); % Peak SLL vs num elements 

AveSLL_neCarray_size)=AveSLL_neCarray_size)+meanCSLLC2;lengthCAl))); % Ave SLL 

vs num elements 

end % array_size 


end % runs 

plotCl;array_size, HPBW_neCl;array_size)/nruns, , l;array_size,theta3db) 

tl=['Random Planar Array ', num2strCnum_elements), ' Elements']; 

t2='HPBW vs Array Size'; 

titleC{tl,t2}); 

xlabelC'Array Size \lambda') 

ylabelC'Half Power Beam Width '^o') 

grid off 

legendC 'Random Array' , 'Analytic') 

%axisCC 1 30 0 360]); 

figure; 

plotCl;array_size, Gain_neCl;array_size)/nruns,.. . 

l;array_size, Gain_neCl;array_size)/nruns - PeakSLL_neCl;array_size)/nruns, ... 
l;array_size, Gain_neCl;array_size)/nruns - AveSLL_neCl;array_size)/nruns, ... 
l;array_size, Gain_neCl;array_size)/nruns-13) % 13 = -101ogl0C20 elements) 
tl=['Random Planar Array ', num2strCnum_elements), ' Elements']; 
t2='Maximum Gain vs Array Size'; 
titleC{tl,t2}); 
axisC [1 20 0 23]); 
grid off 

xlabelC'Array Size \lambda'); 
ylabelC 'Power CdB)'); 

legendC'Max Gain', 'Peak SLL', 'Ave SLL','Ave SLL-Theory' , 'Location' , 'Northwest') 
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C, Random Planar Array: Maximum Gain vs Number of Elements 

This code was used in Figure 15. 

% This code computes far-field gain pattern 
% for planar array composed of N vertical half-wave 
% dipoles at wavelength normalized coordinates [xCn),yCn)] 

% with complex ICn) currents. 

% 

% By M. A. Morgan 29 Oct 07 
% 

% Modified by C.E. Taylor Dec 07 to use Random Planar Array 
% ** This version calculates gain vs number of elements ** 

% 

% This version uses Uniformly current amplitudes CIn=l for all elem) 

% 

clear all 
clf reset 

% Setting up 3-D pattern grid 
Ntheta=91; Nphi=361; 

theta=linspaceC0,pi/2,Ntheta); phi=linspaceC-pi,pi,Nphi); 

[Phi,Theta]=meshgridCphi,theta); 

CT^cosCTheta); ST^sinCTheta); CP^cosCPhi); SP=sinCPhi); 

Rx=ST.*CP; Ry=ST.*SP; Rz=CT; % unit radial vector components 

% Half-wave dipole pattern for Iin=lA (eps handles theta=0 indeterminent) 

E0 = 60*cosCpi*CT./2)./CST+eps); 

Gain_ne = zerosCl,50); % Gain vs number elements array 

PeakSLL_ne = zerosCl,50); % Peak SLL vs num elements 
AveSLL_ne = zerosCl,50); % Ave SLL vs num elements 

% -Build Random Array- 

lambda = 1.0; 
k = 2*pi/lambda; 

array_size = 5.0*lambda; % in wavelengths (lambda) square grid 


nruns = 75; 
for run=l;nruns; 

fprintfC ' run; %2.0f\n', run); 

for num_elements = 5:1:50; 

Nx = num_elements; 

Ny = num_elements; 
sepn_min = 0; 

xp = array_size * randCl,num_elements); 
yp = array_size * rand(l,num_elements); 

% -Stats of Random Array- 

min_spacing = 0.1*lambda; % l/10th wavelength 


while (sepn_min< min_spacing) % min spacing = .1 wavelength 
rp = [xp;yp]'; 

sepn = pdistCrp); % dist between all elements in Random Array 
sepn_min = min(sepn); % find minimum separation 

zsepn = squareform(sepn); % convert linear sepn array to square 

[Rmin Cmin] = find(zsepn==sepn_min,l, 'first'); %returns index 
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fprintfC['Min separation is %6.3f wavelengths,'... 

'between elements %2.0f and %2.0f \n'],... 
sepn_min, Rmin, Cmin) 

if Csepn_min< min_spacing) % then Move one of them 
xpCRmin)=array_size * rand; 
ypCRmin)=array_size * rand; 
end 
end 

fprintfC 'Element %2.0f Located at %8.3f, %8.3f\n',Rmin,rpCRmin,;) ) 
fprintfC 'Element %2.0f Located at %8.3f %8.3f\n',Cmin,rpCCmin,;) ) 

10=1; % assuming uniform current amplitudes 

%phid0=inputC'Enter peak gain azimuth angle in degrees from x-axis: '); 
phid0 = 0; 

phi0=pi*phid0/180; ^convert to radians 

% Phasing ICNx,Ny) array elements to point array in phi0 direction 
Ip= 10 * expC-j*k*CcosCphi0)*xp+sinCphi0)*yp)); 

% Summing far-field pattern 

F=zerosCsizeCE0)); 
for n=l;num_elements 
xpRx = xpCn)*Rx; 
ypRy = ypCn)*Ry; 

F = F + IpCn)* expCj*k*CxpRx+ypRy)); 
end 

% Normalized power pattern 

F2 = absCCE0.*F)/Cnum_elements)).^2; 

DTheta=thetaC2)-thetaCl); DPhi=phiC2)-phiCl); 

IntF2=DTheta*DPhi*sumCsumCCF2).*sinCTheta))); 

Gain=4*pi*F2/IntF2; 

Gmax=maxCmaxCGain)); GmaxdB=10*logl0CGmax); 

% Horizontal pattern (the x-y plane) and 3dB beamwidth 
n=findCphi==0); 
m=findCtheta==pi/2); 

G3dB=Gmax/2; 

GH=GainCm,;); 

mf=findCGH >= G3dB,l, 'first'); 
mp=findCGH >= G3dB,l, 'last'); 

BWH=180*CphiCmp)-phiCmf))/pi; 

GHdB=10*logl0CGH); 


% - calc Side Lobe levels - 

BW = imregionalmaxCGHdB); %MATLAB rtn that finds local maxima->BW array "l"s 

A1 = findCBW); %index pointers to local Maxima 

[bb,ix]=sortCGHdBCfindCBW)), 'descend'); % sort by highest peak->lowest 

SLL= bbCl)-bb; % Side Lobe levels, down from peak 

SLi= AlCix); % index to peaks 

SLa= 180*phiCSLi)/pi; % Angle where peaks occur 

fprintfC'First Sidelobe Level: %6.2f below max at phi= %4.0f \n',SLLC2),SLaC2)) 
fprintfC'Second Sidelobe Level: %6.2f below max at phi= %4.0f \n',SLLC3),SLaC3)) 
fprintfC'Mean Sidelobe Level: %6.2f below max \n',meanCSLLC2:lengthCAl)))) 
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Gain_neCnum_elements)= Gain_neCnum_elements) + GmaxdB; 

PeakSLL_neCnum_elements)= PeakSLL_neCnum_elements)+SLLC2); % Peak SLL vs num 

elements 

AveSLL_neCnum_elements)=AveSLL_neCnum_elements)+meanCSLLC2;lengthCAl))); % Ave 

SLL vs num elements 


end % num_elements 


end % runs 


ne = l:num_elements; 
Gain_ne=Gain_ne/nruns; 

PSSL = Gain_ne-PeakSLL_ne/nruns; % 
ASSL = Gain_ne-AveSLL_ne/nruns; % 
Gth= Gain_ne+ 10*logl0Cl./ne); % 
Gain_th2 = 10*logl0Cne)+10*logl0C1.6); 


Peak SLL actual 

Ave SLL actual 

Theory ASLL (Steinberg) 

% Theory; Gain of N dipoles, incl gain of dipole 


plot(ne, Gain_ne, ne,PSSL,... 

ne,ASSL, ne,Gth, ne, Gain_th2 ) 

tl=['Random Planar Array ', num2str(array_size), ' \lambda square']; 
t2='Maximum Gain vs Number of Elements'; 
titleC{tl,t2}); 
axisC [5 50 0 23]); 

xlabelC 'Number of Elements in Array'); 
ylabelC 'Power (dB)'); 

legendC'Gain' , 'Peak SLL', 'Ave SLL', 'Theory-Ave SLL',... 

'Theory-N dipoles' , ' Location' , 'Northwest' ) 
grid off 


D, DOA determination for the Uniform Linear Array, with two incident signals. 

This code used in Figure 23, as an example of the SDMA technique. 


% Modified by Chris Taylor Jan 2008 to include pn seq generation or Walsh 
% codes, samples per chip, and two incident beams 

% This is a simulation to model the Elam SDMA technique 
% Developed by Frank B. Gross, November 11, 2004 

% define the sample baseband modulation initially as a IkHz phase modulation tone 
% WALSH FUNCTIONS WORK BETTER THAN PN SEQUENCES,Can also use Gold codes 
clear; clc; 

f=lE3; 

T=l/C4*f); 

Nchips=64; 

SPC =4; % Samples Per Chip (4 min) 

tend=SPC*Nchips; % time snapshot end 

t=[0:Ctend-1)]*T/Ctend-1); 
zeta=sin(2*pi*f*t); % modulation m 

%m=zerosCl,lengthCm)); 

% plot array factor for linear array 

N=ll; % number of receive elements to give 10 degree beamwidth 
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d=.5; % spacing in wavelengths between adjacent array elements 


% generate pseudo-random chip codes T long to modulate the 
% received signals from each antenna 

%Construct a PN object usnign MATLAB pn seq gen 
% h = seqgen.pnCShift', 0); 

% % Output chipped PN bits 
% setCh, 'NumBitsOut', Nchips); 

% for i=l:N 

% BCi,;) = generateCh); 

% end 

% B = DmapCB,!); % map 0,1 to +-1 

% use Walsh Sequence instead of pure pn 
H = hadamardCNchips); 
for i=l;N 

BCi,;)=HCi+l , % don't want 1st row 
end 

% spread B by SPC 

B = repmatC B,SPC,1); % repeat by spread factor 

B = reshapeCB,N,SPC*Nchips); % like chips together 
B = B*pi/2; % pi/2 since inside exponential 

%test SDMA concept by elimnating the pn sequence 
%B=zerosCN,SPC*Nchips); 

test=-30*pi/180; 

test2=45*pi/180; 

Vr=zerosC361,lengthCt)); 

Vr2=zerosC361,lengthCt)); 
er=zerosC361,lengthCt)); 
for k=l;361; 

anCk)=Ck-l)-180; 

thCk)=anCk)*pi/180; 

% calculate the actual received signal for each angle using the N element 
% receiver and sum all N signals, 
for i=l;N % create combined signal at angle k 

VrCk,;)= VrCk,;)+ expClj*Czeta+2*pi*Ci“l)*d*sinCtest) +BCi,:))); 
erCk,;)= erCk,;)+ expClj*C 2*pi*Ci“l)*d*sinCthCk))+BCi,:))); 
Vr2Ck,;)=Vr2Ck,;)+expClj*Czeta+2*pi*Ci“l)*d*sinCtest2)+BCi,:))); 


end 

VrCk,:)=VrCk,:)+Vr2Ck,:); 

RlCk)=sumCVrCk,;).* erCk,;)); 

R2Ck)=sumCVrCk,:).*conjCerCk,:))); 
phaselCk,:)=unwrapCangleCRlCI<))); 
phase2Ck,;)=unwrapCangleCR2Ck))); 
end 

AF=sinCN*pi*d*CsinCth)-sinCtest))+eps)./CN*sinCpi*d*CsinCth)-sinCtest)))+eps); 

AF2=sinCN*pi*d*CsinCth)-sinCtest2))+eps)./CN*sinCpi*d*CsinCth)-sinCtest2)))+eps); 

maxAF^maxCabsCAF)); 

maxAF2=maxCabsCAF2)); 

maxR2=maxCabsCR2)); 

maxRl^maxCabsCRl)); 

figure; 
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plotCan,absCR2)/maxR2, 'b' ,an,absCAF+AF2), 'g' ) 

% plotCan,absCR2)/maxR2,'b',an,absCAF),'g',an,absCAF2),'r') 

xlabel('\theta') 

ylabelC'Correlation Magnitude') 

titlestringO = ['Uniform Linear Array - SDMA']; 

titlestring=['N=',num2strCN), ' d=\lambda/2,' , '\theta inc = '... 

,num2strCtest*180/pi), ',' ,num2strCtest2*180/pi) ]; 
titlestring2=['Nchips=' ,num2strCNchips), ' Samples/chip =', num2strCSPC)]; 
titleC{titlestring0,titlestring,titlestring2}) 
axisC[-90 90 0 1]) 
legendC 'Rk' , 'Array Factor' ); 
grid on 

% figure; 

% titlestring=['Correlation Comparison \theta_k =',num2strCtest*180/pi)]; 
% titlestring2='comparison of instantaneous and averaged phases'; 

% plotCt,zeta,t,phase2C121,:)) 

% xlabelC'time (t)') 

% ylabelC'phase') 

% axisC [0 tClengthCt)) -2 2]) 

% grid on 

% titleC{titlestring,titlestring2}); 

% legendC'phase zetaCt)','average phase') 


E, SDMA Technique Applied to the Random Planar Array 

This code was used to generate Figure 24 the Signal DOA determination for a 
random planar array. 

% Chris Taylor Modification to Dr. Gross' SDMA sim 
% to allow for Random Planar array vice ULA 
% Feb 08 Naval Postgraduate School, Monterey CA 
% 

% Algorithm from Frank B. Gross "Smart Antennas for Wireless 
% Communications", c 2005, Chapter 8, p253-260 
% 

% This is a simulation to model the Elam SDMA technique 
% Developed by Frank B. Gross, November 11, 2004 
% WALSH FUNCTIONS WORK BETTER THAN PN SEQUENCES 
% Can also use Gold codes 
clear; clc; clf reset; 

% - Define Constants - 


define the sample baseband modulation initially as a 300MHz phase modulation tone 


f=3E8; 

T=l/C4*f); 

Nchips=64; 

SPC = 16; 

num_elements = 30; 
phi_inc=30*pi/180; 


freq of baseband modulation 
time bandwidth product = 4 
Number of Chips / Time seq 
Samples Per Chip 
number of antenna elements 
incident signal of interest 


lambda = 3E8/f; 
k = 2*pi/lambda; 


% c/f 

% wavenumber 


tend=SPC*Nchips; 

t=[0;Ctend-1)]*T/Ctend-1); 

zeta=sinC2*pi*f*t); 


% time snapshot end 
% discritization of signal 

% modulation zeta (the signal phase sent) 


87 




Ntheta=91; Nphi=361; % Set up angle matrices 

%Ntheta=5; Nphi=13; 

theta = linspaceCO,pi/2,Ntheta); 

phi = linspaceC-pi,pi,Nphi); 

DTheta= thetaC2)-thetaCl); DPhi = phiC2)-phiCl); 

[Phi,Theta] = meshgridCphi,theta); 

CT=cosCTheta); ST=sinCTheta); CP=cosCPhi); SP=sinCPhi); 

Rx=ST.*CP; Ry=ST.*SP; Rz=CT; % unit radial vector components 

% - generate pseudo-random chip codes - 

% T long to modulate the received signals from each antenna 
% "codes" here = "Beta" nomenclature from text 
H = hadamardCNchips); % use Walsh Sequence 
for i=l;num_elements 

codesCi,:)=HCi+l,;); % don't want 1st row 

end 

% spread codes by SPC 

codes = repmatC codes,SPC,1); % repeat by spread factor 

codes = reshapeCcodes,num_elements,SPC*Nchips); % move like chips together 

codes = codes*pi/2; % pi/2 since inside exponential 


Build Random Array 


array_size = 5.0*lambda; % in wavelengths (lambda) square grid 

min_spacing = 0.1*lambda; % l/10th wavelength 

sepn_min =0; % initialize prior to 1st test 

xp = array_size * randCl,num_elements); 

yp = array_size * randCl,num_elements); 

% -Stats of Random Array- 

while (sepn_min< min_spacing) % min spacing = .1 wavelength 

rp = [xp;yp]'; 

sepn = pdistCrp); % dist between all elements in Random Array 
sepn_min = min(sepn); % find minimum separation 

zsepn = squareform(sepn); % convert linear sepn array to square 

[Rmin Cmin] = findCzsepn==sepn_min,l, 'first'); %returns index 
fprintf(['Min separation is %6.3f wavelengths,'... 

'between elements %2.0f and %2.0f \n '],... 
sepn_min, Rmin, Cmin) 

if (sepn_min< min_spacing) % then Move one of them 
xp(Rmin)=array_size * rand; 
yp(Rmin)=array_size * rand; 
end 
end 

fprintfC 'Element %2.0f Located at %8.3f, %8.3f\n' ,Rmin,rp(Rmin,;) ) 
fprintfC 'Element %2.0f Located at %8.3f %8.3f\n' ,Cmin,rp(Cmin,;) ) 

% -Plot Random Array- 

subplotC2,2,l) 

setCgcf, 'PaperOrientation' , 'portrait'); 

plotCxp,yp, 'MarkerSize' ,15, 'Marker' ,, 'LineWidth' ,2,... 

'LineStyle' , 'none' , 'Color' ,[0 10]); 
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tstrg = 'Sensor Cluster, Uniform Excitation'; 
tstrg2= [num2strCnum_elements), ' Elements']; 
titleC{tstrg,tstrg2}, 'FontSize' ,14) 
axisC[0 array_size 0 array_size]) 
axis square 

setCgca, 'XTick' ,linspaceC0,array_size,6)) 
setCgca, 'YTick' ,linspaceC0,array_size,6)) 
xlabelC'x - wavelengths') 
ylabelC'y - wavelengths') 
grid on; 

for i = l;num_elements; 

textCxpCi),ypCi), num2strCi), 'FontSize' ,10) 
end 


% -plot coord system- 

subplotC2,2,2) 

coord_syst_diagCl,phi_inc*180/pi,90,array_size) 
plot3Cxp-array_size/2,yp-array_size/2,zerosClengthCxp)), ... 

'Marker''LineStyle' , 'none' , 'Color' ,[0 1 0]); 
plot3CCPCl,:),SPC1,:), zerosCsizefCPCl,:))), 'c') 
titleC Coordinate System and Scan Plane'); 



% Calc Array Factor for Planar Random Array 

% _ 


% Half-wave dipole pattern for Iin=lA (eps handles theta=0 indeterminent) 
E0 = 60*cosCpi*CT./2)./CST+eps); 

10=1; % assuming uniform current amplitudes 

Ip= 10 * expC-j*k*CcosCphi_inc)*xp+sinCphi_inc)*yp)); 

% Summing far-field pattern 

AF=zerosCsizeCE0)); 
for n=l;num_elements 
xpRx = xpCn)*Rx; 
ypRy = ypCn)*Ry; 

AF = AF + IpCn)* expCj*k*CxpRx+ypRy)); 
end 


%. 

% 

% 

%. 


Calc Received Signals via SDMA 


yr=zerosCNphi,lengthCt)); 
ye=zerosCNphi,lengthft)); 

m=findCtheta==pi/2); % want row for xy plane 

Rxm = Rxfm,;); % IxNphi vector 

Rym = RyCm,;); 


for p=l;Nphi; 

% calculate the actual received signal for each angle using the N element 
% receiver and sum all N signals, 
for n=l;num_elements 
xpRx = xpCn)*RxmCp); 
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ypRy = ypCn)*RymCp); 

re = cosCphi_inc)*xpCn)+sinCphi_inc)*ypCn); 
yrCp,;)=yrCp,;)+expCj*Czeta+ k*CxpRx+ypRy) +codesCn,;))); 
yeCp,:)=yeCp,:)+expCj*C k*Cre) +codesCn,;))); 

end 

RCp)=sumCyrCp,;).*conjCyeCp,:))); 
end 

maxR=maxCabsCR)); 


Calc Gain, 


Side Lobe Levels, HPBW, etc. 


% Horizontal pattern (the x-y plane) and 3dB beamwidth 
n=findCphi==0); 
m=findCtheta==pi/2); 

F2 = absCCE0.*AF)/Cnum_elements)).^2; % Normalized power pattern 
IntF2=DTheta*DPhi*sumCsumCCF2).*sinCTheta))); 

Gain=4*pi*F2/IntF2; 

Gmax=maxCmaxCGain)); GmaxdB=10*logl0CGmax); 

% Convert Correlation Magnitude to Gain Pattern 
RF2 = absCCE0Cm,;).*R)/Cnum_elements)).^2; 

GainR2 = 4*pi*RF2/IntF2; %use same normalizing factor 
GRdB = 10*logl0CGainR2); 

% Steinberg's Random element Gain pattern Eqn 8.10 
F2_Rand_ave = F2.*Cl-l/num_elements)+l/num_elements; 

Gain_Rand_ave = 4*pi*F2_Rand_ave/IntF2; 
GRA_dB=10*logl0CGain_Rand_aveCm,:)); 


G3dB=Gmax/2; 

GH^GainCm,;); 

mf^findCGH >= G3dB,l, 'first'); 
mp=findCGH >= G3dB,l, 'last'); 

BWH=180*CphiCmp)-phiCmf))/pi; 

GHdB=10*logl0CGH); 

AFmax = maxCabsCAFCm,;))); 

% - cqIc; Side Lobe levels - 

BW = imregionalmaxCGHdB); %MATLAB rtn that finds local maxima->BW array "l"s 

A1 = findCBW); %index pointers to local Maxima 

[bb,ix]=sortCGHdBCfindCBW)), 'descend'); % sort by highest peak->lowest 

SLL= bbCl)-bb; % Side Lobe levels, down from peak 

SLi= AlCix); % index to peaks 

SLa= 180*phiCSLi)/pi; % Angle where peaks occur 

fprintfC First Sidelobe Level: %6.2f below max at phi= %4.0f \n' ,SLLC2),SLaC2)) 
fprintfC Second Sidelobe Level: %6.2f below max at phi= %4.0f \n' ,SLLC3),SLaC3)) 
fprintfC 'Mean Sidelobe Level: %6.2f below max \n' ,meanCSLLC2:lengthCAl)))) 


% _ 

% 

% Plot Array Factor vs. SDMA Pattern 

% _ 
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% subplotC2,2,3:4) 

% plotCphi*180/pi,absCR)/maxR,'b',phi*180/pi,absCAFCm,;))/AFmax,'g') 

% xlabelC'Xphi') 

% ylabelCCorrelation Magnitude') 

% tstring=['Number Elements=',num2strCnum_elements),... 

% ' Nchips=',num2strCNchips)]; 

% tstring2=['\phi inc = ',num2strCphi_inc*180/pi), ' M=', num2strClengthCs))]; 

% titleC{tstring,tstring2},'FontSize',14) 

% axisC[-180 180 0 1]) 

% legendC'Rk','Array Factor'); 

% grid on 


%. 


% Plot Gain vs. SDMA Gain Pattern 

% _ 


GHMax=maxCGHdB); 

GRMax^maxCGRdB); 

GRAMax=maxCGRA_dB); 

Diff=GRMax-GHMax; % SDMA only calculated in one planar cut vice sphere,... 

% Normalization factor is adjusted to acct for diff. 

phid=phi*180/pi; 
subplotC2,2,3:4) 

plotCphid,GRdB-Diff, 'b' ,phid,GHdB, 'g' , 'LineWidth' ,1.5) 

%plotCphid,GHdB,'b',phid,GRdB-Diff,'g',phid,GRA_dB,'r','LineWidth',2) 

v=axis; vCl)=-180; vC2)=180; if vC3) < -60, vC3)=-60; end 

axisCv) 

tl = ['Horizontal Plane: G_{max} = ' ,num2strCGmaxdB, '%8.1f'),'dB'] ; 
t2 = ['HPBW= ', num2strCBWH, '%8.1f'), ' deg',... 

' SLL=' , num2strCSLLC2), '%6.2f below Peak')]; 
titleC{tl,t2}, 'FontSize' ,14) 
xlabelC'Xphi (deg)' , 'FontSize' ,14) 
ylabelC 'Gain Pattern CdB)' , 'FontSize' ,14) 
legendC'SDMA Correlation', 'Planar Array Gain'); 
grid on 
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F, Incident Signal Phase, Walsh-Hadamard Codes & Received Signal 

This code used to generate Figure 22 showing the received signal phase, Walsh- 
Hadamard Spreading Codes and the real part of the chipped signal r(t). 

% This is a simulation to model the Elam SDMA technique 
% Developed by Frank B. Gross, November 11, 2004 

% define the sample baseband modulation initially as a IkHz phase modulation tone 
% WALSH FUNCTIONS WORK BETTER THAN PN SEQUENCES,Can also use Gold codes 

% Modified by Chris Taylor Jan 2008 to include pn seq generation or Walsh 
% codes, samples per chip, and two incident beams 

% This short code shows the phase of incident signal zetaCt), Walsh codes 
% and the Real part of the received signal r_nCt) 

clear; clc; 

f=lE3; 

T=l/C4*f); 

Nchips=64; 

SPC =4; % Samples Per Chip (4 min) 

tend=SPC*Nchips; % time snapshot end 

t=[0;Ctend-1)]*T/Ctend-1); 
zeta=sinC2*pi*f*t); % modulation m 

% plot array factor for linear array 

N=ll; % number of receive elements to give 10 degree beamwidth 
d=.5; % spacing in wavelengths between adjacent array elements 

% generate pseudo-random chip codes T long to modulate the 
% received signals from each antenna 

% use Walsh Sequence instead of pure pn 
H = hadamardCNchips); 
for i=l;N 

BCi,;)=HCi+l,;); % don't want 1st row 

end 

% spread B by SPC 

B = repmatC B,SPC,1); % repeat by spread factor 

B = reshapeCB,N,SPC*Nchips); % like chips together 
B = B*pi/2; % pi/2 since inside exponential 

figure; 

subplotC3,l,l) 
plotCt,zeta, 'b') 
axisC [0 T -2 2]) 
xlabelC'time (t)') 
ylabelC 'Phase - \zetaCt)') 
grid on 

subplotC3,l,2) 

plotCt,BCll,;), 'r') % use arbitrary WH row 

axisC [0 T -2 2]) 

xlabelC'time (t)') 

ylabelC 'Walsh-Hadamard Code') 
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grid on 


subplotC3,l,3) 

plotCt,realCexpCj*Czeta+BCll,:)))), 'm') 

axisC [0 T -2 2]) 

xlabelC'time (t)') 

ylabelC'Real Part of r_nCt)') 

grid on 
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